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APPROXIMATE MEASUREMENT OF COLOR TEMPERATURE 
WITH A BARRIER-LAYER PHOTOCELL! 


By W. E. KNOwLES MIDDLETON AND A. R. RAMSEY? 


Abstract 
An objective method of measuring color temperature is described, using 
apparatus found in nearly every photometric laboratory. 


Introduction 


The measurement of color temperature by direct visual comparison with 
a standard is one of the more exacting photometric measurements. It is 
not surprising, therefore, that photoelectric methods have been devised 
(6, 10) that are capable of very high accuracy if carefully applied. Their 
general principle involves simultaneous or successive measurements of the 
radiation from the test source by two receivers of different spectral sensitivity, 
as for example two photocells with different cathode surfaces, or one photocell 
with two interchangeable filters. 

No novelty in principle is claimed for the present note, the purpose of which 
is to direct attention to the fact that objective measurements of color tem- 


perature can be made with equipment that is in regular use in many photo- 
metric laboratories, and to report an investigation into the errors of such 
measurements. 


Equipment 

On the North American continent, it is becoming almost standard practice 
to perform routine photometry by means of a barrier-layer cell and a so-called 
‘zero-resistance’ circuit (1, 7). The cell is, of course, fitted with a filter that 
corrects its spectral sensitivity to that of the normal eye (2) with a sufficient 
degree of accuracy for the purpose in hand. The zero-resistance circuit 
gives adequate linearity for most purposes, and, if a cell with a reasonably 
low drift can be selected, a precision nearly equal to that of visual photo- 
metry can be attained, with an immense saving of time. The circuit used 
in the present work is essentially the same as that of Projector, Laufer, and 
Douglas (7), with the addition of a galvanometer reversing key along the 
lines described by Wenner (9), which effectively doubles the sensitivity of 
the galvanometer. This key (K, in Fig. 1), incidentally, consists of two 
S.P.D.T. ‘Micro-switches’ mounted side by side and operated by a common 
button, giving a very quick make and break. No thermoelectric effects due 
to this key have been detected with the portable galvanometer used. 


1 Manuscript received October 10, 1947. 
Contribution from the Division of Physics and Electrical Engineering, National Research 
Laboratories, Ottawa, Canada. Issued as N.R.C. No. 1705. 
2 Present address: Northern Electric Company, Montreal, Que. 
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The slide wire is of the Kohlrausch type, as generally used in such cir- 
cuits, and can be read to the nearest 1/10,000 of its length with great ease. 
The entire equipment, in fact, is equivalent to that found in many industrial 
photometric laboratories, with the addition of a suitably diaphragmed holder 
for the barrier-layer cell, and two filters. It is best to use the entire sensitive 
surface of the cell, using the diaphragms merely to screen off stray light. 





- |+ 
STD. CELL 





SLIDE WIRE 





SENSITIVITY DAMPING 
Fic. 1. ‘Zero-resistance’ circuit used. 


Cell and Filters 

The cell used in our experiments is The General Electric Co. light-sensitive 
cell, supplied with a green filter by the makers. The combination has a spectral 
sensitivity very similar to the normal I.C.I. luminosity curve, though small 
variations occur from cell to cell. In Fig. 2 are plotted the relative spectral 
sensitivity of the cell used (full line) and the normal luminosity curve (dashed 
line). For the present purpose an exact correspondence is not at all necessary. 

The spectral transmittance of the pair of red and blue filters used was 
measured on an automatic recording spectrophotometer and is shown in 
Fig. 3. These filters were chosen from among a large number of available 
types, on the following general principles. 
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Fic. 2. Spectral sensitivity of barrier-layer cell with green filter (full line) compared with 
the spectral luminosity function for the ‘normal’ I.C.I. observer (dashed line). 








Fic. 3. Spectral transmittance of red filter (Corning 3480, 5.36 mm. thick) and of blue 
filter (Corning 4308, 5.11 mm. thick). 
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(1) They should divide the entire visible spectrum fairly evenly between 


them. 
(2) Their transmittance should have as regular a slope as possible. 
No doubt many other pairs of filters would have been equally suitable. 


Calibration 


Because of the difficulty of measuring accurately the spectral sensitivity of 
the photocell, it is better to calibrate the apparatus by the use of a color 
temperature substandard, a lamp for which the relation between voltage 
and color temperature is known. Such lamps are usually calibrated in the 
range 2000° to 2850° K. Light approximately equivalent to sources of higher 
temperature can be obtained by the use of blue filters (3, 4) with the color 
temperature substandard. 

With the lamp and the photocell in positions that will give a convenient 
illumination, say about 20 I. ft.-?, readings are taken with the filters in position 
in the order red, blue, red, blue, red. The mean of the two readings with the 
blue filter J, , is divided by the mean of the three readings with the red, J,, 
and the resulting ratio, R, is adopted as an index of the color temperature of 
the source. This procedure is repeated with the source set to give light of 
various color temperatures within the desired range. The result of such a 
calibration is shown in Fig. 4. The points obtained by means of a blue glass 








o. 
2000 2400 2800 3200 3600 °K 4000 


Fic. 4. Relation between color temperature T and ratio R. The points indicated by squares 
refer to ‘color temperatures’ obtained by means of a blue photometric filter. 
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filter (NBS 15430) lie on a curve nearly parallel to the curve obtained for 
the unfiltered tungsten lamp. If it is decided to adopt the scale of color 
temperatures furnished by the combination of lamp and filter, there is no 
reason why the method should not be used if it is calibrated with the particular 
combination employed.* 

Users of the ‘zero-resistance’ circuit will recognize that there is a limited 
range of illuminations within which it is both linear and sensitive. In taking 
the measurements at low values of color temperature, when little illumination 
is available from the substandard, a certain amount of sensitivity will have 
to be sacrificed because of the comparatively low slide wire readings; but 
above about 2250° K., where readings of several hundred divisions can easily 
be obtained, the precision of a single measurement of R, made as above, is 
at least +0.002, corresponding to about +8°K. This assumes that the 
lamp voltage is regulated sufficiently well, as is usual in photometric labora- 
tories. It is, however, advisable to operate within a 3 to 1 range of illumina- 
tion if such precision is desired; if a tolerance of 20° K. is permissible, a range 
of 10 to 1 in illumination can be used. But since the apparatus also serves 
to measure the illumination, and the lamp is usually mounted on a movable 
carriage, it is just as easy to adopt a standard value of illumination and work 
near it. 

The measurement is not very sensitive to changes in the temperature of 
the filters. “Experiments in a cold chamber at about 0° C. showed an over- 
estimation of about 20° K. at 2850° K., when a calibration made at 25° C. 
was used. This overestimation decreased to almost zero at 2500° K. It 
may be concluded that the effect is less than one degree in color temperature 
per degree Centigrade, and thus unimportant in the ordinary laboratory, 
unless the filters are exposed to strong illumination and warmed up in this 
way. 

For routine checks it is adequate to take only three readings, red, blue, and 
red. With only ordinary care a measurement of this sort should be reliable 
to within +15° K. if the illumination does not differ by a factor of more than 


3 from that used in.the calibration. 


Limitations of the Method 


The most obvious limitation of the method is that it demands that the 
illumination be kept constant while three readings (or possibly only two) are 
made. This limits the method to the photometric laboratory, for which it 

‘was obviously designed. A recent German patent (8) describes a method 
of overcoming this difficulty by the use of two barrier-layer photocells, but 
some experiments by the authors lead them to doubt whether this method 
would be sufficiently accurate for laboratory purposes. 


* The ‘color temperatures’ corresponding to these points were obtained by the method of Judd (5). 
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The second and not so obvious limitation of the method is that it will 
assign a color temperature to any illumination that is presented to it, even if 
this is of a color to which it is highly inappropriate to apply the idea of color 
temperature. This limitation applies to all objective methods. The user 
will have to decide for himself, for example, whether the color temperature 
of fluorescent lamps should be measured in this way. The fundamental 
method, direct visual comparison in a photometric field, makes it very evident 
whether a given illumination lies on the Planckian locus. 

The authors have not found the ‘drift’ of the barrier-layer cell at all serious 
when working with an illumination of about 20 to 301. ft.-? Any drift that 
occurs after the first minute or so has a negligible effect on the color temperature 
measurement. The situation is different where an illumination of 200 I. ft.-? 
is in question; here the effect of drift is at least as great as all the other errors 
together, and such high values of illumination should not be employed in these 
measurements. 
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PERMEABILITY STUDIES 
III. SURFACE AREA MEASUREMENTS OF CARBON BLACKS! 


By J. C. ARNELL AND G. O. HENNEBERRY 


Abstract 


The modified Kozeny equation has been found to be satisfactory for the 
measurement of the specific surfaces of carbon blacks having average particle 
diameters ranging from 0.01 to 0.1 w to within +10%. Comparative data 
were obtained from electron microscope counting and from low temperature 
nitrogen adsorption isotherms. The three methods examined gave results that 
were in satisfactory agreement, except when the carbon black was porous, and 
then the adsorption value was extremely large. 


Introduction 


In two preceding papers (2, 3), a method has been described for the measure- 
ment of the specific surface of fine powders by permeability methods. This 
method depends on the use of a modified Kozeny equation which may be 
written as: . 








G= shear | €Pug gP 43x 3 ee | (1) 
LS,(1 — €) | konS,(1 — €) 3\V 7M 
where G = _ rate of gas flow in cubic centimeters per second, 
A = cross-sectional area of the bed in square centimeters, 
h = pressure head in centimeters of mercury, 
€ = volume of pore space per unit volume (porosity), 
@ = fraction effective void area, 
Pi = p/p = density of the gas at unit pressure, 
L = _ length of the bed in centimeters, 
S, = specific surface of the powder in square centimeters per cubic 
centimeter, 
Pug= density of mercury in grams per cubic centimeter, 
g = acceleration due to gravity, 
P = (p; + p2)/2 = mean pressure in centimeters of mercury, 
ky = _ shape factor, approximately equal to 5a/e, 
n = viscosity of the gas in poises, 
6 = variable factor, having a value of approximately 0.9, 
Ro = gas constant, 
T = absolute temperature, and 
M = molecular weight of the gas. 


This equation has been tested against the data obtained from the flow of 
air through beds of finely ground quartz and inorganic pigment powders and 
found to give satisfactory results. 


1 Manuscript received September 5, 1947. 
Contribution from the Defence Research Chemical Laboratories, Ottawa, Canada, 
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In this paper, the method has been extended to a series of carbon black 
samples, and the resulting specific surfaces have been compared with those 
obtained from electron microscopy and low temperature nitrogen adsorption. 


Apparatus and Technique 


The first permeability measurements on the carbon blacks were made with 
the apparatus that has been described previously (3), but the low rate of 
flow of air through the packed beds rendered this apparatus impractical and 
a new apparatus was designed. 

A schematic drawing of the new apparatus is shown in Fig. 1. The Lucite 
tube containing the packed bed of powder, C, was connected into the apparatus 
by means of rubber pressure tubing, and a pressure differential across the bed 





Fic. 1. Diagram of permeability apparatus. 


was maintained by means of a falling bead of mercury in the capillary tube, 
B, of uniform bore. The short manometer, A, provided a means of measuring 
the pressure head across the bed, while the large manometer, E, was used to 
measure the absolute pressure in the apparatus. The five liter bulb, D, was 
introduced as a surge volume to compensate for the fact that a larger volume 
of gas left the bed than entered it, and, without a large surge volume, the flow 
of gas through the bed resulted in an absolute pressure increase in the 
apparatus. 

The technique of measuring permeabilities was as follows. With the packed 
tube in position, stopcock S_: was closed and stopcock S; was opened to vacuum. 
As the system was evacuated, a pressure differential was built up across the 
bed and finally the mercury bead was drawn up the capillary tube. When the 





omen 
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mercury bead was near the top of the capillary tube, stopcock S; was closed 
and the measurement of the rate of flow could be made. The mercury bead 
was allowed to fall about 7 cm. to establish equilibrium conditions, and then 
the time of fall of the bead for each centimeter was recorded over a distance 
of at least 20cm. In general, the elapsed time taken by the mercury bead to 
fall 20 cm. was used in the flow calculations. On completion of a run, stop- 
cock S: was opened and the pressures in the apparatus equalized. If the same 
operation were repeated, a second run could be made at a lower mean pressure. 

From a knowledge of the cross-sectional area of the capillary tube, B, and 
the pressures in the system, the mass flow at the fixed pressure head could be 
calculated by means of the equation: 


AopiPid 
C= oe, (2) 
t 
where Ao = cross-sectional area of the capillary tube in centimeters, 
P,; = pressure on high pressure side of bed in centimeters of mercury, 
d = _ length of fall of mercury bead in centimeters, and 


t = time necessary for bead to fall d cm., in seconds. 

It was found necessary to follow the fall of the mercury bead over a con- 
siderable distance, because the movement of the bead was somewhat erratic. 
This was thought to be due to surface forces between the mercury and the 
glass, even though the glass had been thoroughly cleaned in hot acid and the 
mercury acid-washed and distilled before the apparatus was assembled. 
This nonuniformity of fall of the mercury bead appeared to be generally 
nullified over a distance of 10 cm., and the time of fall over two consecutive 
10 cm. lengths usually agreed to within several seconds. 

In the apparatus used for these experiments, the capillary tube had a cross- 
sectional area of 0.0547 sq. cm. and the mercury bead was 5.0 cm. long. 


Theoretical 


With all powders previously studied, the relation between the mass flow and 
the mean pressure could be expressed by the equation: 


A 


G=5 (BP+C), (3) 
where B = slope of the line resulting from a G vs. P plot, and 
C = corresponding intercept. ' 


In the paper immediately preceding this one, it was pointed out that the 
specific surface of a powder could be calculated from the intercept value by 
using only the second half of Equation (1), or by writing :— 


s, = Ee ee. (4) 
3'N 7M C(1 — €) 


where 6 is taken as 0.9. 


The majority of the carbon blacks studied had average particle sizes that 
were about an order of magnitude smaller than those of the finest powders 
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previously studied, and, when the mass flow through a carbon black was 
measured at different mean pressures, it was generally found that the mass 
flow remained constant to within the experimental error. At first it was 
thought that this was an indication that pure Knudsen flow was occurring 
through the beds, but subsequently several beds of larger blacks with high 
porosities were found to give flows of the type expressed by Equation (3) with 
very small values for the slope. Calculations of the Adzumi pore radii (1) 
were made for those beds that had a finite slope, when plotted according to 
Equation (3), and the values of the radii fell in the region of 2-5 X 10-* cm. 

The Knudsen semiempirical equation for the flow of gas through a capillary 
tube (8, p. 296) was examined for the flow through a capillary having a radius 
of 5 X 10-* cm. and a length of 1 cm. This equation takes the form: 


_ & eer Asa ee (5) 
hp 8mL '1+aP3N M L’ 





T 


where c; = 2.00 ’ Ro éeee 2 AT : R, and R = radius of the capillary in 


centimeters. 

The resultant curve is shown in Fig. 2, and it is interesting to note that the 
minimum occurs at a pressure of just over 30 cm. of mercury and that in the 
pressure range of the experiments described here (10 to 65 cm. of mercury), 
the total variation is only 2%. It therefore seemed reasonable to assume 
that the region of the flow of air through the beds of carbon black was near 
the minimum in the Knudsen curve and to take an average value through 
this region for calculating the values of the specific surface. This average 
value was found to be equivalent to giving 6 a value of 0.96 in Equation (4). 

All the values of the specific surface to be given in the experimental section of 
this paper have been calculated by means of Equation (4) and where the G vs. P 
plot of the data for any bed was found to have a slope, 6 was taken as 0.9 and 
when the mass flow, G, appeared to be constant, 6 was given the value of 0.96. 


1.175 





sec. X 10% 


L, C286 SP. 





P, cm. mercury 
Fic. 2. Plot of Knudsen equation for a fine capillary. 
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Experimental Results 


A series of 12 commercial carbon blacks were chosen as representative of 
the different methods of manufacture. The specific surface was determined 
using at least three beds of each black, and the mass flow was generally 
measured at seven mean pressures ranging from 15 cm. of mercury to 70 cm. 
of mercury. 

The runs were all made at room temperature and the porosities of the beds 
were calculated using densities calculated by X-ray diffraction. 

In order to check the method, the surface areas of the carbon blacks were 
also calculated from particle counts made from electron microscope photo- 
graphs by means of the equation: 


So= nd? 


; (6) 


where K is a shape factor, which was taken as 6. 1, since all the blacks appeared 
spherical (for perfect spheres, K = 6.0). Low temperature nitrogen adsorp- 
tion isotherms were also determined and the surface areas were then calculated 
by means of the Brunauer, Emmett, and Teller equation (4) and also by the 
Harkins and Jura equation (6). 

The experimental results are all listed in Table I, together with pertinent 
data pertaining to each black. 

















TABLE I 
SURFACE AREAS OF CARBON BLACKS BY PERMEABILITY, ELECTRON MICROSCOPE, AND ADSORPTION 
METHODS 
Surface area, m.?/cc. 
Black Type* | Porosity N2 adsorption |Butane** 











Permeability | Electron —| adsorp- 
microscope} BET | H & J tion 


Thermax MT -- o 13:6 14.5 15.8 19 

Kosmos 20 SRF 0.554 43) 
0.568 43> 42 47 41 52 29 

0.579 41 

Kosmos 40 HMF 0.584 67) 
0.634 59> 65 72 78 82 45 

0.638 69} 

Shawinigan CF 0.729 67 | 
100% Comp. 0.762 82, 74 78 94 — 85 

0.791 72) 

Statex B FT 0.601 97) 

; 0.610 99| 
0.619 103+ 99 72 99 87 80 

0.636 98 | 

0.645 96) 

Statex A CF 0.603 98 | 
0.639 100+ 100 86 136 152 78 


0.656 | 103) 
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TABLE I—Concluded 


SURFACE AREAS OF CARBON BLACKS BY PERMEABILITY, ELECTRON MICROSCOPE, AND ADSORPTION 
METHODS—Concluded 





Surface area, m.?/cc. 


Black Type* | Porosity N2 adsorption |Butane** 
Permeability | Electron |—————-————| adsorp- 


microscope] BET |'H & J tion 


Kosmobile 77 EPC 0.638 155) 


0.670 | 150] 
0.672 | 155| 146 161 193 202 80 
0.673 | 131/ 
0.674 | 131! 
0.675 | 155} 
0.715 93 
0.721 70 

Dixie densed MPC 0.696 164) 

HM 0.704 | 144} 150 165 228 269 94 

0.716 | 141) 

Kosmobile S HPC 0.653 178) 
0.666 | 197} 182 162 277 304 105 
0.673 171) 
0.692 | 126 
0.707 71 
0.718 45 

Spheron 9 EPC 0.676 189) 
0.696 | 180} 184 144 223 225 80 
0.696 | 182} 

Voltex cl 0.707 344) 
0.717 | 328} 326 220 =e 1320 281 
0.765 | 306! 

Neo Spectra Ink 0.823 | 435) 

Mk. II Black 0.825 | 418$ 428 259 1410 1960 = 

0.828 | 430) 











* This refers to the type of process by which each black was made, and further details can be 
found in the manufacturer's literature. 7 

** Data taken from Table 18, page 90, Developments and Status of Carbon Black by Isaac 
Drogie. United Carbon Company, Inc., Maryland, U.S.A., 1945. 


Discussion 


A comparison between these independent methods of measuring the specific 
surfaces of a number of carbon blacks is provided by Table I. It is of interest 
to examine these methods critically and to consider the relative merits of each. 


Electron Microscope 

From a cursory consideration, the electron microscope appears to provide 
the most direct method of evaluating the specific surface of a fine powder. 
There are, however, a number of indeterminate factors which reduce the 
accuracy of the method considerably and which also apply, in general, to 
similar measurements with larger particles using the optical microscope. 





SOR ea Parte MTS 


FRSA: 
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The method, which is based on Equation (6), is a statistical one and con- 
sequently requires that the diameters of a large number of particles be 
measured. This is particularly important if there is a large spread between 
the diameters of the smallest and the largest particles. In general, a group 
of between 5 and 10 electron microscope photographs will have enough 
sufficiently dispersed particles to permit the measurement of several hundred 
diameters. These are not always sufficient data to ensure a good statistical 
average, as the following example indicates. In a preceding paper (2), a 
sample of finely ground quartz powder was studied and the diameters of 335 
quartz particles were measured on electron microscope photographs. Two 
of these particles were found to be at least twice as large as any other particle 
in the group and about 30 times as large as the particles in the most probable 
size range. The specific surface c: this sample of quartz was found to be 
6 X 10' sq. cm. per cc. and, if the two large particles were discarded, the 
value of the specific surface was approximately doubled. While this is an 
extreme case, it seems to illustrate the care necessary in using microscopic 
methods for powders having a wide range of particle sizes. 


The shape factor in Equation (6) is another indeterminate source of error, 
as it involves the surface-volume ratio of the particles. As long as the 
particles are nearly spherical, as was the case with the carbon blacks being 
considered here, a value can be assigned for the shape factor, which has no 
greater an error than 2 to 3%. As the shape of the particles diverges from 
that of a sphere, the error in the assigned value of the shape factor increases 
until it may reach 20 to 30% in the case of needlelike particles. 


With many powders it is very difficult to get a good dispersion of the parti- 
cles, and agglomerates are always present in the field of vision. When the 
powder being examined has a wide range of particle sizes, it is often very diffi- 
cult to decide between an individual large particle and an agglomeration of 
small particles. In the case of the carbon blacks, it was felt that some of the 
large particles may have been discarded for this reason and the values of the 
specific surface, in consequence, may be slightly high. 

Another source of error, that applies in the main to particles of the size of 
the finest carbon blacks, is the difficulty of enlarging the particles to a size that 
permits fairly accurate counting. The carbon blacks were all enlarged photo- 
graphically at around 17,000 diameters and the negatives were then examined 
under a visual enlarger, which raised the magnification to about 125,000. 
Even at this power, a fine black was composed of particles with an average 
apparent diameter of about 2 mm. For example, the Neo Spectra Mk. II 
carbon black had the following particle distribution: 


Lessthan0.01 0.01-0.02 0.02-0.03 0.03-0.04u Total 


3 65 64 5 137 
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These particles were too small to be subdivided into any more classes and it is 
apparent that slight errors in the measuring scale, which was based on the 
final viewing magnification, could markedly change the distribution, and, in 


consequence, the specific surface. 


Adsorption 

The adsorption method measures any internal surface that is open to 
the adsorbate, in addition to the external surface of the particles. This 
results in a value of the surface area of a porous material higher than the values 
of the specific surface determined by methods based on the determination of 
external surface alone, as is the case with the microscopic and the permeability 
methods. There is no way of telling a priori when a powder consists of porous 
particles and consequently, if only the external surface is desired, the adsorp- 
tion method is not reliable. The results in Table I suggest that the majority 
of the carbon blacks studied are either slightly porous or have roughened 
surfaces, and the last two blacks in the table are quite obviously very porous, 
with large internal surfaces. 

There are two theories that may be applied to a low temperature nitrogen 
isotherm in order to calculate the surface area of the adsorbent. The first 
of these is the Brunauer-Emmett-Teller (BET) method (4), which is valid 
for a sigmoid isotherm and provides a means of estimating the volume adsorbed 
in a monolayer and of calculating the surface area, if an area is assigned to 
an adsorbed molecule. The method probably has an error of at least 10%, 
and, as Gregg (5) has pointed out in an excellent discussion of the method, the 
error may run as high as 50% if the isotherm is not a sigmoid (BET Type I1) 
isotherm. The surface area may also be calculated from the semiempirical 
isotherm equation of Harkins and Jura (6), which appears to fit a greater 
number of isotherm shapes than the BET equation. With the carbon blacks 
there was one case (Shawinigan) where the Harkins and Jura plot was too 
curved for a calculation to be made and the isotherms of the Voltex*and the 
Neo Spectra resembled the Langmuir type, and consequently the BET 
equation failed to give a straight line plot, but, in general, the specific surfaces 
of the carbon blacks calculated using the two equations were in fair agreement. 
Zettlemoyer and Walker (11) also report satisfactory agreement for the sur- 
face areas of a number of samples of active magnesia calculated using the two 
equations. When the two equations lead to surface area values that show 
reasonable agreement, the average value of the surface probably is correct 
within 10%. 

The butane adsorption values of the surface area listed in Table I are 
included only as a matter of interest. These values are average figures that 
have been quoted in the manufacturer’s literature and are only intended to 
give the order of magnitude of the surface. When it is considered that 
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surface areas calculated from butane adsorption isotherms are generally con- 
siderably smaller than those calculated from nitrogen adsorption isotherms, 
the carbon black samples used for this study appear, in general, to have been 
fairly representative. 


Permeability 

A modification of the Kozeny equation for the calculation of the surface 
area of a powder from gas-permeability measurements has been proposed by 
Rigden (9, 10) and discussed further by Lea and Nurse (7). This modification 
is similar to that proposed by one of the authors (2, 3) except for one significant 
point. Rigden has used the assumption that the fraction void area is equal 
to the fraction void volume ¢. This is theoretically true. However, in the 
case of packed powders, where the dimensions of the individual channels 
through the bed are varying continually, it has been considered unlikely that 
the effective void area would approach the theoretical one and, in consequence, 
a term defined as the fraction effective cross-sectional area, a, has been used 
instead of the theoretical fraction void area, €, in the derivation of Equation (1). 
The use of a has been discussed at some length in the preceding papers of this 
series and an empirical relation between @ and € has been used for all specific 
surface calculations to date. This relation is: 


log a = 1.41e — 1.40. (7) 


In general, Equation (1) and the Rigden equation should give values of the 
specific surface that lie within the experimental error of the method. 

There are several indeterminate factors that introduce errors into the 
permeability method. One of these is the difficulty of ensuring that a bed is 
homogeneous and uniformly packed. Compaction of the bed in increments 
appears to be the recommended method. However, Lea and Nurse (7) have 
shown that beds compacted in one stroke give specific surfaces that are not 
more than 3% lower than an increment-packed bed. This suggests that a 
well packed bed is reasonably homogeneous. 

Lea and Nurse also point out that there is a general decrease in the specific 
surface values with increasing porosity using the Rigden equation. The use 
of Equation (1), involving a, appears to eliminate this fall-off over the inter- 
mediate porosity range, although it is still present at the higher porosities. 
This fall-off at the high porosities may be partially due to the fact that high 
porosity beds are less likely to be homogeneous and the particles tend to bridge, 
forming an occasional large channel. Such a channel would act as a major 
conduit and prevent many of the smaller channels from playing their part, 
with a resultant apparent decrease in calculated specific surface. Rigden (10) 
also reports a fall-off in the value of the specific surface, calculated for the same 
bed, with decreasing mean pressure, which is an effect that has never been 
observed by the authors. 
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With soft materials there is always the possibility of deforming the particles 
and creating appreciable areas of contact between particles. As this deform- 
ation would be a function of compacting pressure and would appear as a 
decrease in specific surface, it may be assumed to be absent as long as the 
calculated specific surface does not increase with increasing porosity. 


Conclusions 


It would appear that each of the methods examined is probably subject 
to an error of about 10%, when applied to such small particles. Other 
methods, such as sedimentation, could not be applied owing to the fineness of 
the carbon blacks. The permeability method is considered to be the most 
satisfactory of the methods examined, as it provides a relatively quick estimate 
of the external surface of very fine powders, and in the final analysis has about 
the same accuracy as the other methods, which are more time-consuming and 
require fairly elaborate apparatus. 
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THE EFFECT OF VARIATIONS IN THE GENERAL COMPOSITION 
OF SAMPLES IN SPECTROGRAPHIC ANALYSIS! 


By G. O. LANGSTROTH AND D. ANDRYCHUK? 


Abstract 


Variations that occur in intensity ratios belonging to pairs of elements on the 
addition of various extraneous substances to standard samples have been studied 
using a common type of condensed spark discharge. The extent of the variations 
depended on the test elements chosen, the nature of the added substance, the 
line pair investigated, and on whether or not the standard samples were buffered. 
The variations bore no apparent relation to the relative ionization potentials or 
atomic weights of the test elements, nor to any physical property of the added 
substances. In general, they were smaller for the line pairs having the most 
nearly equal excitation potentials. Usually, but not always, they were smaller 
for buffered than for unbuffered samples. Illumination of the electrodes with 
ultraviolet light, though resulting in steadier operation of the discharge, or 
alterations of the capacity in the sparking circuit did not reduce the variations 
materially. Operation of the discharge at constant current coupled with the use 
of a buffer however did reduce the variations considerably for the instances 
studied. 

It has been commonly recognized that changes in the general composition 
of samples introduced into an arc or spark source may cause variations in 
the relative intensities with which spectral lines of two elements under test 
are radiated, even though the elements are always present in constant amounts. 
This effect has an important bearing on the accuracy of spectrographic analyses 
of samples varving considerably in composition. A study (4) of intensity 
ratio changes in the radiation from a condensed spark when various substances 
were added to samples containing magnesium, lead, and cadmium provided 
no obvious clue to the factors responsible for the variations. An extension 
of this study, involving a much greater range of data, provides the basis for 
the present article. 

Apparatus and Procedure 

The condensed spark discharge was produced by a circuit consisting of a 
0.014 uf. condenser, a 30 wh. inductance, a thermocouple ammeter, and a 
spark gap in series. The condenser was connected across the secondary of a 
transformer (full rating, 20 kv.) whose primary voltage was controlled by an 
autotransformer. The electrodes, of the point and movable plane type (1), 
were copper. The plane electrode on which the sample was dried was 1.2 cm. 
by 1.0cm. in area, and an electrode separation of 3.1 mm. was used throughout 
the work. 

Spectra were photographed with a medium quartz spectrograph, using a 
slit width of 45 uw. An image of the discharge was formed on the collimating 
lens by a quartz condensing lens at the slit, so that light from all regions of 
the discharge contributed to the intensity of spectral lines. A standard 

1 Manuscript received August 26, 1947. 


Contribution from the Department of Physics, University of Alberta, Edmonton, Alta. 
(The experimental work was done in the Physics Department, University of Manitoba.) 


2 Holder of a Bursary under the National Research Council of Canada at the University of 
Manitoba. 
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exposure of 60 sec. was used; this was made feasible by means of a moving 
wire screen weakener (2) before the source. For one sample containing a 
complex salt, determination of intensity ratios of lines of the two anion con- 
stituents was made with an antimony weakener (3) covering about half the 
spectrograph slit, in order that strong lines in the weakened section of the 
photograph might be comparable in density to fainter lines in the unweakened 
section. 

The plane electrode was loaded by delivering 0.050 cc. of the desired liquid 
sample to it with a special pipette calibrated to 2%, and drying under low 
pressure. As a general practice, more electrodes were prepared than were 
needed, and those with uniform dried deposits were selected for sparking. 
During an exposure the plane electrode was moved so that the spark traversed . 
a zigzag pattern over it, continually striking to a fresh surface of the deposit. 
For a large part of the work the spark was operated at the minimum breakdown 
voltage, the primary voltage being adjusted until the discharge was just on 
the point of stopping; for some special experiments to be mentioned later 
other conditions were employed. 


The standard samples whose general composition was varied by the addition 
of extraneous substances consisted of eight solutions each containing three 
test elements (Table I). Constituents (e.g., silver in Sample 1) that would 

TABLE I 


COMPOSITION OF THE STANDARD SAMPLES 








Salt Salt Salt 
Sample aie eine 
rs Conc., ms Conc., - Conc., 
Kind mgm./cc. Kind mgm./cc. Kind mgm./cc. 

1 AgNO; 0.194 | SnCl. . 6H.0 0.94 MgCl, . 6H,0 1.05 
2 AgNO; 0.194 | SnCl. . 6H20 0.94 PbNO; : 2.00 
3 MnCl, . 4H,O0 1.50 NiCl . 6H,O 1.00 MgCl. . 6H,O 1.05 
4 NiCl . 6H.O 1.00 | SnCl. .6H,O 0.94 MgCl. . 6H:O0 1.05 
5 BaCl, . 2H2O 8.00 | SnClh . 6H2,O 0.94 CdCl; . 2H,O 4.00 
6 SrCl. . 6H,O 12.0 AICI; . 6H,0 1.00 MgCl, . 6H,0 1.05 
1.60 


7 SrCl. . 6H:0 10.0 NiCl . 6H,0 3.00 | ZnCl, 
8 KAI(SOx)2 . 12H2O 4.00 MgSO,.7H,O} 4.00 - 


have been precipitated if made up with the other two were prepared in separate 
solutions and loaded on the electrodes separately. The test elements were 
chosen with consideration for ionization potential, atomic weight, and posses- 
sion of suitable spectral lines (Tables III and IV). The concentrations 
required to yield comparable intensities for the selected lines were found by 
trial. 
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Each of the standard samples 1 to 7 was divided into five portions, of which 
one was left unaltered and each of the others had a different extraneous 
material added. The added materials, their concentrations, and the minimum 
ionization potentials associated with them, are given in Table II. Sample 8 
was treated in a similar way but the added materials were different. Duplicate 


TABLE II 


EXTRANEOUS SUBSTANCES ADDED TO THE STANDARD SAMPLES 








To Samples 1 to 7 To Sample 8 
. Conc., | Min. atomic ion. = : Conc., | Min. atomic ion. 

Substance mgm./cc. pot., ev. Substance mgm./cc. pot., ev. 
KCl 6.00 4.32 (K) CdSO, 12.0 8.96 (Cd) 
NH.Cl 1.00 13.0 (Cl) (NH4)2SO4 30.0 10.3 (S) 
LiCl 3.60 5.36 (Li) Na2SO, 15.0 15.12 (Na) 
H;BO3 10.3 8.26 (B) ZnSO, 15.0 9.36 (Zn) 

TABLE III 


IONIZATION POTENTIAL, V;, AND ATOMIC WEIGHT, M, OF THE TEST ELEMENTS 





Element Vi, ev. M Element Vi, ev. M 
K 4.32 39 Mn 7.41 55 
Ba 5.19 137 Ag 7.54 108 
Sr 5.67 87 || Mg 7.61 24 
Al 5.96 27 Ni 7.61 59 
Sn 7.30 118 Cd 8.96 112 
Pb 7.38 207 Zn 9.36 65 





spectra for all portions of a standard sample were photographed on a single 
plate, along with a calibration spectrum taken with a rotating step-sector and 
a brass arc operated under constant conditions. A knowledge of the relative 
intensities of lines of the calibration spectrum was unnecessary, since informa- 
tion on the changes in intensity ratios only was required and each plate con- 
tained a complete set of data fora sample. The intensities of selected lines of 
the test elements (in arbitrary units, not necessarily the same for each) were 
determined in the usual way with a microphotometer (5). Typical line 
ratios studied are given in Table IV along with data on their relative excitation 
potentials. 


The experiments outlined above were repeated in exactly the same way 
with the added feature that each of the standard samples was ‘buffered’ with 
50 mgm. per cc. of potassium acetate. For the unbuffered samples the load 
of salts on the electrode ranged from 9.1 X 10-7 gm. per mm.? for Sample 1 
with no addition, to 104 X 10-7 gm. per mm.? for Sample 7 with boric acid 
added; for the buffered samples it varied from 217 X 10-7 to 269 X 10-7 gm. 
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per mm. for the same instances. In these experiments, the amount of 
extraneous material added was not varied. Only in one instance (addition 
of boric acid to unbuffered Sample 4) was any attempt made to investigate 
the effect of changing the amount of the added extraneous substance. 


It was noted during the experiments outlined above that, when the spark 
was operated at minimum breakdown voltage, additions of extraneous material 
caused changes in the spark circuit current. A subsidiary experiment was 


TABLE IV 


TYPICAL LINE PAIRS, THE SYMBOLS USED TO REFER TO THEM, AND THE RATIO OF THEIR 
EXCITATION POTENTIALS, V/V’ 








. Symbol A Symbol B Symbol C Symbol D 
Element eee EN et sh 6 en Pa = Ss 
Pair A/r’ | V/V" | A/N’ V/V’ A/V’ V/V’ A/N’ V/V’ 
Ag/Sn | 3383/3262 | 0.75 3281/3262 | 0.78 | 3383/3175 | 0.85 | 3281/3175 | 0.87 
Ag/Pb | 3383/2614 | 0.64 3281/2614 | 0.66 | 3383/2833 | 0.84 | 3281/2833 | 0.86 
Sn/Pb 3175/2614 | 0.75 | 3262/2614 | 0.85 | 3175/2833 | 0.99 | 3175/3640 | 0.99 


Sn/Mg 3175/3838 | 0.73 3262/3838 | 0.82 | 3175/2852 | 1.00 | 3262/2852 | 1.12 
Ni/Sn 3525/3262 | 0.73 | 3515/3262 | 0.75 | 3525/3175 | 0.82 | 3415/3175 | 0.85 
Ni/Mg 3515/3838 | 0.61 3415/3838 | 0.62 | 3515/2852 | 0.83 | 3415/2852 | 0.84 
Ni/Mn | 3415/2610 | 3515/2610 3515/4033*| 1.18 | 3415/4033*] 1.19 
Mn/Mg | 4033*/2852 | 0.71 2632/2852 — — — — 
Sr/Al 4832/3092 | 1.09 4812/3082 | 1.10 | 4742/3961 | 1.40 | 4722/3944 | 1.40 


| 

Ag/Mg | 3383 2a | 0.5 | 3281/3838 | 0.64 | 3383/2852 | 0.84 | 3281/2852 | 0.87 
| 
| 


Sr/Mg 4832/3838 | 0.74 | 4742/3838 | 0.74 | 4722/2852 | 1.01 | 4812/2852 | 1.02 
Al/Mg 3961/3838 | 0.53 3082/3838 | 0.68 | 3944/2852 | 0.72 | 3092/2852 | 0.93 











Ni/Sr | 3525/4742 | 0.81 3515/4812 | 0.84 | 3415/4742 | 0.84 | 3493/4832 | 0.85 
Ni/Zn | 3525/3345 | 0.46 3515/3345 | 0.47 | 3415/3345 | 0.47 | 3493/3345 | 0.47 
Sr/Zn 4832/3345 | 0.56 | 4742/3345 | 0.57 | — — - — 
Cd/Sn 4679/2840 | 1.33 | 4679/3175 | 1.48 | 3404/3175 | 1.70 | 3404/2863 | 1.70 
Ba/Sn 3501/2840 | 0.74 | 4903/3175 | 0.97 | 4903/2863 | 0.97 | 4823/3175 | 1.00 
Cd/Ba | 4679/4283 | 1.47 4679/4903 | 1.51 | 4679/3501 | 1.78 | 3403/3501 | 2.09 
K/Al 4047/3082 | 0.76 4047/3961 | 0.98 _— — — 











* A blend of 4033 and 4034 A. 


therefore performed with Sample 4 plus lithium chloride, in which intensity 
ratios were determined at different currents as obtained by adjustment of 
the transformer primary voltage and read from the thermocouple ammeter. 
Variations with current were observed for some ratios but not for others. 
For example, with a buffered sample the nickel-magnesium ratios C and D 
(Table IV) decreased by about 16% on increasing the current from 1.3 to 
1.8 amp., and by about 25% on increasing it to 2.2 amp.; on the other hand the 
ratios A and B remained effectively unchanged. With an unbuffered sample 
the former ratios decreased by about 14% on increasing the current from 1.5 
to 2.0 amp. and the latter ratios remained essentially constant. Accordingly, 
the main experiments with Samples 1 and 2, with and without the buffer, 
were repeated, the spark being operated at constant current rather than 


minimum breakdown voltage. 
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It was further noted during the experiments that some of the added materials 
caused changes in the minimum breakdown voltage. Accordingly, the main 
experiments on unbuffered Sample 1 were repeated with the additional feature 
that the plane electrode immediately under the point was illuminated with 
intense ultraviolet light (the image of an iron arc, run at 4.0 amp. and 34 cm. 
distant, was formed on the electrode by an f : 4 quartz combination). This 
arrangement reduced the breakdown voltage and resulted in a steadier opera- 
tion of the spark. <A further complete set of data for Sample 1 was obtained 
under this condition with the capacity in the spark circuit increased from 0.014 


to 0.026 uf. 
. B Results 


The results for the main experiments with the spark operated at minimum 
breakdown voltage are given in Fig. 1 for Samples 1 to 7. The symbols A, 
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Fic. 1. Variations in intensity ratios caused by addition of extraneous materials to the 
samples, as obtained with the spark operated at minimum breakdown voltage. The manner of 
presentation 1s described in the text. 








44 CANADIAN JOURNAL OF RESEARCH. VOL. 26, SEC. A. 


B, C, and D refer to the various line pairs (Table IV), which are given in order 
of increasing excitation potential ratio from left to right for each pair of ele- 
ments. The vertical lines, whose length represents the intensity ratio obtained 
on addition of an extraneous substance divided by that obtained for the 
unaltered sample, are plotted in groups of four for each line pair; the members 
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Fic. 1 (continued). Variations in intensity ratios caused by addition of ‘extraneous 
materials to the samples, as obtained with the spark operated at minimum breakdown voltage. 
The manner of presentation is described in the text. 


of each group from left to right refer respectively to additions of potassium 
chloride, ammonium chloride, lithium chloride, and boric acid. The solid 
dots denote the corresponding values obtained when the samples were 
buffered with potassium acetate (they are referred to the intensity ratio 
for the unaltered buffered sample). If the addition of an extraneous substance 
had no effect on intensity ratios the ends of the vertical lines and the dots 
would have ordinate values of unity. The plotted values are averages of 
nine sets of data for Sample 1, four sets for Samples 2, 3, and 4, and two sets 
for Samples 5, 6, and 7. Results for the same element pair obtained with 
different samples are given separately (e.g., silver-tin in Samples 1 and 2) 
since the composition of the load on the electrode differed. 
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r The effectiveness of various amounts of added boric acid in changing the 
> intensity ratios for unbuffered Sample 4 is indicated in Fig. 2. In general, 
d boric acid was found to be more effective than the other added substances in 
e ; changing intensity ratios. 
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Fic. 2. Changes in intensity ratios caused by the addition of various amounts of boric acid 
to unbuffered Sample 4. The symbols refer to the following line pairs: Ni/Sn, O 3525/2863, 
+3415/3262; Ni/Mg, O 3515/3838 (A), + 3415/2852 (D); Sn/Mg, @ 3175/3838 (A), 
+ 3175/2852 (C), O 3262/2852 (D). 


Average results obtained for unbuffered Sample 8 with the spark operated 
at minimum breakdown voltage are shown in Fig. 3. The method of presenta- 
| tion is similar to that of Fig. 1. The numbers of each group of four solid lines 


K /Al 
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° Fic. 3. Variations in intensity ratios on the addition of extraneous material to unbuffered 

) Sample 8, which contained a complex salt (potassium aluminium sulphate.) 

refer from left to right to values of the intensity ratios (relative to that for 
. the unaltered sample) obtained on addition of zinc sulphate, sodium sulphate, 


cadmium sulphate, and ammonium sulphate respectively. The dashed lines 
refer to the values for additions of potassium chloride, ammonium chloride, 
) lithium chloride, and boric acid to unbuffered Sample 6, and have been trans- 
ferred from Fig. 1 for comparison. 












46 CANADIAN JOURNAL OF RESEARCH. VOL. 26, SEC. A. 


A comparison of the average results obtained with various modifications of 
the sparking procedure is given in Fig. 4. The method of presentation 
is similar to that of Fig. 1, with results obtained by each modification of the 
procedure denoted by a special type of line. 
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Fic. 4. A comparison of the variations in intensity ratios caused by addition of extraneous 
materials, as obtained by various modifications of the procedure. The solid lines refer to values 
obtained with the spark operated at minimum breakdown voltage and have been transferred from 
Fig. 1 for comparison; the dotted lines denote values obtained with constant spark current, 
the broken lines those obtained with ultraviolet illumination on the electrode, and the dashed 
lines those obtained with increased capacity and ultraviolet illumination. 


Discussion 


Examination of Fig. 1 shows that extremely large variations may be caused 
in intensity ratios by changes in the general composition of samples. It is 
evident that the magnitude of the variation depended on the line pair and 
element pair under consideration, on the nature of the substance added to the 
sample, and on whether or not the sample contained a buffer substance. 
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Furthermore, it appears probable from the results given in Fig. 2 that for 
some line pairs the size of the variation depends strongly on the amount of the 
extraneous substance added, but for others it does not. These findings are 
in agreement with earlier results (4). 


The data of Fig. 3 indicate that large intensity ratio variations are caused 
by added substances even when the test elements of a pair are in combination 
in a complex salt and so at least must have the same distribution over the 
electrode surface and be equally subject to the impinging spark discharge. 
The aluminium/magnesium data of Fig. 3 show that additions of zinc, sodium, 
cadmium, or ammonium sulphates produced variations comparable to those 
caused by additions of potassium, ammonium, or lithium chlorides or boric 
acid. 

The data of Fig. 4 indicate that illumination of the electrode with ultra- 
violet light, although producing a steadier discharge, did little to reduce 
silver/tin or silver/magnesium intensity ratio variations caused by the 
addition of extraneous material; this was true also when the capacity of the 
spark circuit was approximately doubled. 


Operation of the spark at constant current rather than minimum breakdown 
voltage appears to offer some advantage, since with buffered samples it led to 
considerably reduced intensity ratio variations for silver/tin, silver/lead, and 
silver/magnesium (Fig. 4). With unbuffered samples, however, the variations 
were similar. to those obtained under minimum breakdown voltage operation. 


The remainder of the discussion deals with an analysis of the data of 
the main experiments on Samples 1 to 7 with minimum breakdown voltage 
operation. The element pairs nickel/tin, silver/lead, tin/lead, silver/tin, 
nickel/manganese, strontium/aluminium, nickel/magnesium, silver/magne- 
sium and tin/magnesium form a series having nearly constant ionization 
potential ratios (ranging from 0.96 to 1.04), but atomic weight ratios graded 
from 0.50 to 4.90. The average variations in intensity ratio, considering 
line pairs A and D and all added extraneous substances, were compared for 
the members of this series in the attempt to find some trend involving the 
atomic weight ratio. No significant trend was evident. For example, for 
an atomic weight ratio of 0.50 (nickel/tin) the average variation with unbuf- 
fered samples was 8%, but for a ratio of 0.52 (silver/lead) it was 51%; for 
ratios of 0.91 (silver/tin) and 4.9 (tin/magnesium) the average variations 
were 22% and 31% respectively. 

A similar comparison was made of the average variations for the series 
aluminium/magnesium, nickel/zinc, nickel/manganese, silver/tin and cad- 
mium/tin. This series has a roughly constant atomic weight ratio (ranging 
from 0.91 to 1.08) but ionization potential ratios graded from 0.70 to 1.23. 
No trend relating the size of the average variation and the ionization potential 
ratio was evident. For example, with a ratio of 0.70 (aluminium/magnesium) 
the average variation was 14%, with a ratio of 1.04 (silver/tin) it was 21%, 
and with a ratio of 1.23 (cadmium/tin) it was 14% for unbuffered samples. 


48 CANADIAN JOURNAL OF RESEARCH. VOL. 26, SEC. A. 


Furthermore, element pairs having ionization potential and atomic weight 
ratios close to unity were associated with average variations as large as those 
for pairs having both ratios quite different from unity. The average variations 
for silver/tin, barium/tin, and cadmium/barium were 22, 22, and 29% 
respectively with unbuffered samples. 

The effectiveness of added substances in producing intensity ratio variations 
could not be related to any of their physical properties. Phe average varia- 
tion, considering line pairs, A, B, C, and D and all element pairs in unbuffered 
Samples 1 to 7, was progressively greater for added substances in the order 
ammonium chloride, potassium chloride, lithium chloride, and boric acid. 
Some element pairs provided exceptions. For example, with nickel/tin the 
order was ammonium chloride, boric acid, lithium chloride, and potassium 
chloride, while with nickel/manganese it was boric acid, ammonium chloride, 
potassium chloride, and lithium chloride. On the average the order with 
buffered samples was potassium chloride, ammonium chloride, lithium 
chloride, and boric acid, but some element pairs provided exceptions. With 
silver/lead the order was boric acid, lithium chloride, potassium chloride, and 
ammonium chloride, and with barium/tin it was lithium chloride, boric acid, 
potassium chloride, and ammonium chloride. 

In general, the intensity ratio variations caused by addition of extraneous 
material were less for buffered than for unbuffered samples. On the average, 
they were about 40% less with potassium chloride as the addition, 42% less 
with lithium chloride, and 57% with boric acid, but with ammonium chloride 
additions there was little difference. Some element pairs provided exceptions, 
e.g., the variations for silver/tin were greater with buffered samples. It was 
noted, however, that reproducibility was higher with buffered samples. 


The average variation for line pairs having the most widely different excita- 
tion potentials was from 40 to 50% greater than that for pairs having the most 
nearly equal excitation potentials, all element pairs and additions being taken 
into account. Considering the various additions separately, it was 63% 
greater on the addition of boric acid, 53% greater on addition .of lithium 
chloride, and 26% greater on addition of ammonium chloride to unbuffered 
samples; and 58, 100, and 17% greater on the addition of these substances to 
buffered samples. It was about the same for both types of line pairs with 
additions of potassium chloride. There were very few individual instances 
in which smaller variations were found for the line pairs having the most 
widely different excitation potentials. They were found for strontium/ 
nickel on the addition of potassium chloride, for strontium/nickel and alumi- 
nium/magnesium on addition of lithium chloride, for aluminium/magnesium 
on addition of boric acid, and for zinc/nickel on addition of ammonium 
chloride, all with unbuffered samples. 

In the light of these results it seems important that application of spectro- 
graphic methods to the analysis of samples that may vary considerably in 
general composition be preceded by a somewhat extensive investigation of the 
variability of determinations with the type of material under consideration. 
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The data provide no clue for the choice of a suitable internal standard element 
since the variability of intensity ratios does not seem to depend in any system- 
atic way on the ionization potential or atomic weight ratio of an element pair. 
There is evidence that in general some advantage is to be gained by choosing 
lines with nearly the same excitation potential, but this feature in itself does 
not ensure stability (see the results for line pairs C of tin/magnesium, C of 
strontium/magnesium, and D of tin/lead, which have excitation potential 
ratios of 1.00, 1.01, and 0.99 respectively). For the most part, the use of a 
buffer substance appears desirable, but in certain instances it may lead to less 
rather than more stability. In some instances at least the use of a buffer and 
operation of the spark at constant current offers considerable advantage. 
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THE THERMAL CONDUCTIVITY OF NAPALM*-GASOLINE GELS! 


By G. O. LANGSTROTH AND F. ZEILER 


Abstract 


A slightly modified form of the rapid and simple method suggested by Hutchin- 
son for the measurement of the thermal conductivity of liquids has been found to 
suffer from convection effects with samples of low viscosity, except with con- 
ductivity tubes of very small diameter. For more viscous liquids such as glycerol 
it was found to be adequate under all conditions studied. The method has been 
applied in the determination of the conductivity of Napalm-gasoline gels. For 
temperatures T between —50° and 50° C., and Napalm concentrations C between 
0 and 10%, the conductivity k in cal. sec. cm.~! per degree C. is described to 
better than 1% by the relation, k X 10° = 29.7 — 0.068 7 + 0.11 C. The 
temperature coefficient of resistance of the unaged tungsten filaments used in the 
tubes differed considerably from the value given by the International Critical 
Tables for aged tungsten filaments. For temperatures 7) between —50° and 
50° C., the coefficient a per degree C. is given to better than 1% by the relation, 
a X 10° = 427 — 1.89 J) + 0.0088 73 — 0.000021 T%. 


A method for measuring the thermal conductivity of liquids, similar to 
that used commonly in the past for measurements on gases, has been sug- 
gested by Hutchinson (6), who has given some data in support of its reliability. 
The method has marked advantages in speed and simplicity over the more 
accurate procedure developed by Bates and his coworkers (1, 2, 3, 4). This 
article describes an application of the method in the determination of the 
conductivity of Napalm-gasoline gels over a wide range of temperature and 


Napalm concentration. In addition, evidence is presented that the method 
as slightly modified by us yields spuriously high values with liquids of low 
viscosity because of the occurrence of convection, except with conductivity 


tubes of small diameter. 
Apparatus and Procedure 


The conductivity tubes (Fig. 1) differed from that described by Hutchinson 
(6) in the use of a more rugged filament (diameter of the spiral 0.020 as 
compared to 0.008 cm.), and the position of one of the side arms. These 
modifications were required by the necessity for filling the tubes with a jelly- 
like substance. The spiral filament, F, was taken from an unused 40 w. 
incandescent lamp and had a resistance of about 25 ohms at 25°C. It was 
spot-welded to No. 18 tungsten leads which were sealed into selected Pyrex 
tubing with the filament under slight tension and as nearly as possible coin- 
cident with the tube axis. External wires were coated with a resistant 
lacquer. For investigating the validity of the method a series of nine tubes 

* Napalm is a commercial product consisting of a basic aluminum soap of naphthenic, oleic, 
and palmitic acids in the approximate ratio 1: 1:2 


1 Manuscript received December 4, 1947. 
Contribution from the Physics Department, University of Alberta, Edmonton, Alberta. 
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with internal diameters ranging from 1.270 to 0.330 cm. was used. The 
length of the filament spiral was from 4.138 to 4.238 cm., and the other 
dimensions were as indicated in the scale drawing of Fig. 1. The external 
diameter of the Pyrex tube for the series ranged from 1.505 to 0.546 cm. 
Tubes were filled to a level about 1 cm. above the upper filament weld, either 
by suction or with the aid of a small syringe. Care was necessary with the 
more viscous gel samples to avoid the inclusion of air bubbles or breakage of 
the filament. 


Fic. 1. 


The tube, with both arms stoppered, was inserted vertically in a constant 
temperature bath. An 8 liter water bath with a special mercury thermostat 
control was used for work above 0° C. Temperature variations were less than 
0.001° C. as indicated by series of readings with a Beckmann differential 
thermometer. Bath temperature was read from a 0.1° mercury thermometer 
that had been calibrated against a certified standard. For work below 0° C., 
a covered, carefully insulated 3 liter alcohol bath was used. A heating coil 
(dissipation about 15 w.), activated by a toluene thermostat control, worked 
against a cooling deviee. This device consisted of a 6.5 cm. diameter metal 
cylinder containing a 4.5 cm. diameter cylinder, both filled with solid carbon 
dioxide in acetone; the cooling effect was controlled by varying the depth of 
insertion of either in the bath. Temperature variations were less than 
8.01° C. as indicated by a 0.1°C. pentane thermometer used to measure 
the bath temperature. No evidence of temperature variation was detected 
with either bath in the electrical measurements described below. 


The circuit used to determine the resistance, R, of the tube filament, the 
current, 7, through it, and the potential drop, V, across it, is represented in 
Fig. 2. The filament, F, and the calibrated box, D, served essentially as 
end coils to the 1 m. slide wire, AB, in a bridge circuit having boxes, S and Q, 
as tatio arms. A mercury rocker switch consisting of two rows of eight cups 
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linked in pairs (shown just above AB in Fig. 2) permitted interchange of 
F and D in the bridge. On interchange the balance point on the slide wire 
was shifted by an amount / cm. such that, 


(R — Ra) = pl, (1) 


where R, denotes the known resistance of D, and p represents the resistance 
per centimeter of the slide wire.* The potential drop, V, was measured with 





to 
potentiometer 


See Se A switch links 
es ae a position 1 


pes ed ee! l position 2 


Fic. 2. 


a Leeds and Northrup Student potentiometer with a certified standard cell 
when the switch was set to Position 2; that across D (Vz) was measured when 
the switch was set to Position 1. The current, 7, was calculated from i = 
V2/Ra; it was adjusted by means of the resistance box, C. 


No variation in the value of p (0.005615 + 0.000003 ohm per cm.) was 
detected for the range of currents used in these experiments. The slide wire 
was calibrated for uniformity along its length. The sensitivity of the galvano- 
meter G (1.3 wV per mm.) was sufficient to permit a balance to the nearest 
0.1 mm. on the slide wire with the minimum filament current. Potential 
drops were measured to 0.05% or better. Filament currents ranged from 
about 0.01 to about 0.05 amp. In general, the bridge was operated with 


500 ohms in S and Q and about 25 ohms in D. 
* Strictly speaking, Equation (1) should contain the resistance difference between the leads to 
F and to D from the switch, but this did not exceed 0.00005 ohm. 
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In meking a determination, a steady current of about 0.01 amp. was passed 
through the filament with the tube inserted in the bath at the desired tem- 
perature. Thermal equilibrium was attained in two to three minutes as 
shown by the constancy of the filament resistance with time. At least four 
readings were made of the shift / in the balance point on interchanging F and 
D, and V and Vg were determined. The procedure was repeated with three 
additional filament currents, requiring in all from 20 to 25 min. The mean pl 
values were plotted against the filament power dissipation, P = V. Va/R., 
and the conductivity, k, was calculated from the slope and the P = 0 intercept 
of the curve in accordance with Equation (7). The magnitudes of the various 
experimental quantities are illustrated by the following figures for an inter- 
mediate filament current. 


Ra = 24.217 + 0.0012 t = 0.02958 + 0.00002 amp. 
Va = 0.7165 + 0.0003 v. pl = 0.3705 + 0.00032 

V = 0.7255 + 0.0003 v. P = 0.02146 + 0.00002 w. 

L = 54.75 + 0.02cm. 


Theoretical Considerations 


If, for the conductivity tube described, the absence of convection and of 
temperature discontinuities at solid—liquid boundaries is assumed, and if it 
is further supposed that conduction of heat from the filament is entirely radial, 
it may be shown easily that for thermal equilibrium 


P = kG(T —T,) = J(Ti — To), (2) 


where P and k are expressed respectively in watts and cal. sec.—'cm.— per degree 
C., and T, T;, and Ty denote the temperatures of the filament, the sample—Pyrex 
interface, and the bath,in °C. The tube constants G and J are described by 


G = 26.30 h/In(r;/r) | 
(3) 


Se 0.0686 h/\n(ro/7r;) 


where 7 and /: denote respectively the radius and length of the filament spiral, in 
cm., and r;and ro represent the internal and external radii of the Pyrex tube. 
In calculating the numerical factor for the latter equation, the conductivity 
of Pyrex was taken to be 0.00261 cal. sec.—'cm.~! per degree C. at 25°C. (7, 
Vol. 5, p. 217); as will appear later, this value need not be known with high 


accuracy. 


Since in practice the temperature of the filament, T, exceeds that of the 
bath, To, by less than 5° C., one may write to a high degree of approximation 


(T — To) = (R = Ro)/aoRo , (4) 
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where R and Rp denote the filament resistance at T and Ty , and a represents 
the temperature coefficient of resistance appropriate to the range 7) to J. It 
follows from Equations (2) and (4) that 


R— Ro 4 


k = P/G(~— = - = 


5 
aoRo J ( ) 


or one may write 
kG kG 
P 


aoRo (1 + kG/J) oe ~ ato(l es kG/J) (6) 


On this basis, a plot of P vs. R should yield a straight line having a P = 0 
intercept of Ro and a slope o such that 


k = @aoRo/G(1 — caoRo/J) (7) 


The terms P/J in Equation (5), RG/J in Equation (6), and gapRo/J in 
Equation (7) are introduced because of the temperature drop across the 
Pyrex walls of the tube, and may be considered as small correction terms. 
J need not be known with high accuracy, but errors in G cause similar errors 
ink. Since r; and r enter Equation (3) logarithmically, errors introduced in 
G by inaccuracies in aligning the filament with the tube axis are of minor 
importance except with very small diameter tubes. The uncertainty in h is 
of more importance. The filament spiral is opened out for a short distance 
near each weld, and it is difficult to judge with high accuracy the value to be 
used for the effective spiral length. For this reason mainly, the error in G 
as calculated from tube geometry may be as high as 4%. Once the validity 
of the basic assumptions has been established, however, more accurate values 
of G may be obtained by standardizing the tubes with suitable liquids of known 
conductivity. 


The Validity of the Method 


The results of some preliminary tests with water suggested the presence of 
convection in the conductivity tubes used. A series of experiments was there- 
fore performed to seek for any possible dependence of single determinations 
on tube diameter, filament power input, and the viscosity of the test substance. 
A value for the ‘apparent’ conductivity was calculated from Equation (5) 
for each determined filament resistance and power input. Since G and J 
were calculated from Equation (3) for this part of the work, the apparent 
conductivities obtained with a given tube may be subject to a systematic 
error of about 4% introduced in specifying tube geometry. Some typical 
data are given in Table I. 

It is evident from Table I that the apparent conductivity of liquids with the 
viscosity of ethyl alcohol or less increased progressively with increasing power 
input when large diameter tubes were used. Even with the smallest power 
input the values obtained were commonly significantly greater than the liter- 
ature values. Furthermore, reduction of the tube diameter caused both of 
these characteristics to appear less strongly; with a diameter of 0.33 cm. the 
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dependence on power input was undetectable and the apparent conductivities 
agreed with accepted values. The conductivities for viscous liquids such as 
glycerol were independent of power input and tube diameter, and agreed with 
the accepted values. The behavior described by the data is typical of effects 
produced by convection, which apparently was present with low viscosity 
liquids in the larger tubes but not with the more viscous samples. The reason 
for the difference between this finding and Hutchinson’s statement (6) that 
the method does not suffer from convection is not clear. 
TABLE I 
THE ‘APPARENT’ CONDUCTIVITY AT 25°C. FOR VARIOUS SUBSTANCES LISTED IN ORDER OF 


INCREASING VISCOSITY 


App. k(X 10°, cal. sec. cm.~! per degree C.) 
Tube for stated power inputs (X 10%, w.; 
Substance diam., expressed approximately) 

cm. 





22 | 38 





Acetone* . 266 57 61 65 
.667 51 49 


Gasoline .016 38 44 46 
.667 37 39 41 
.330 29 28 


Distilled waterT .270 144 150 154 161 
- .330 138 132 136 





Ethyl! alcohol** .270 53 51 56 
.667 47 47 49 
0.330 44 45 45 





Glycerol*t . 266 69 68 68 
.667 70 69 70 
.330 64 68 67 


Napalm-gasoline gel (3.5%) .270 27 28 28 
667 | 28 29 29 











* Laboratory grade; k = 0.00043 c.g.s. units (7, Vol. 5, p. 227). 
tk = 0.00143 c.g.s. units (2). 
** Laboratory grade, denatured; k = 0.00045 and 0.00042 c.g.s. units for water contents of 
5 and 0% respectively (3). 
*+ U. S. P. grade, determined water content less than 2%; k = 0.00069 and 0.00068 c.g.s. 
units for water contents of 2 and 0% respectively (2). 


Aside from the occurrence of convection under certain conditions, the basic 
assumptions stated in the section on theoretical considerations appear to be 
good approximations. Heat losses by conduction through the filament leads 
were less than 0.4% of the power input as estimated on the basis of Gregory 
and Archer's analysis (5). The filament resistance vs. power input curves 
are straight lines (Fig. 3, A) as required by Equation (6), except under con- 
ditions in which convection is serious (Fig. 3, B). Finally, under appropriate 
conditions the method yields results in agreement with those obtained by 
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others using different methods. For example, the average values for the 
conductivity of glycerol at 25° C. calculated from the slopes and intercepts 
of P vs. pl curves obtained with tubes of 1.270, 1.266, 1.248, 1.161, 1.016, 
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Fic. 3. Experimental curves for (A) glycerol and (B) acetone, at 25°C. The number 
attached to each curve refers to the tube used; the internal diameters of tubes, 2, 4, and 5 were 
respectively 1.266, 1.016,and 0.667 cm. Incurve B the symbols © and + represent measure- 
ments taken in the order of increasing and decreasing power input respectively; the broken line 
ts straight. 


0.667, 0.330, and 0.330 cm. internal diameter were respectively 0.000680, 
0.000688, 0.000696, 0.000705, 0.000708, 0.000717, 0.000681, and 0.000665 
cal. sec.'cm.—! per degree C. The conductivity lies between 0.00068 and 
0.00069 (footnote to Table 1). Since a G value found from Equation (3) 
was used in the calculations the above values are subject to the systematic 
error mentioned previously. Standardization of the tubes consisted in 
adjusting the G values so that the data for each tube yielded a conductivity 
of 0.000680 cal. sec.-'cm.—! per degree C. The standardized G values were 
used in the work with Napalm-gasoline gels. 

Application of the theory of the section on theoretical considerations to 
experimental data permits calculation of various quantities not directly 
measured. The following figures for glycerol at 25° C. indicate the magnitude 
of the temperature difference across the sample and across the Pyrex tube 
walls. With a tube of 1.266 cm. diameter, (JT — 7;) and (JT; — To) were 
respectively 0.14° and 0.0014° C. at a power input of 0.002430 w., and 2.12° 
and 0.022°C. at an input of 0.03770 w. The corresponding figures for a tube 
of 0.330 cm. diameter were 0.11° and 0.0045° C. at an input of 0.002544 w. 
and 1.52° and 0.070° C. at an input of 0.03985 w. The tube diameters as 
well as the power inputs are at opposite ends of the range investigated. 


Results for Napalm-gasoline Gels 


The gels were prepared by mixing Napalm with a commercial grade of 
gasoline and maintaining mechanical stirring for half an hour. A gel having 
5% Napalm concentration contained 5 gm. of Napalm per 100 gm. of gel. 
Measurements were made on samples that had been freshly mixed and on 
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samples aged for several months, but since no difference in conductivity was 
detected no distinction has been made in presenting the results. The measure- 
ments on gasoline alone were made with a tube of 0.330 cm. diameter; the 
conductivities agree fairly well with those obtained by extrapolating the k vs. 
Napalm concentration curves to zero concentration (Fig. 4). 


per 


degree C. 


k X 105, cal. sec.~! cm.7 





Napalm concentration, % by weight 


Fic. 4. The thermal conductivity of Napalm-—gasoline gels. The plotted points are expert- 
mental values, the lines were calculated from the empirical relation given in the text. The 
number attached to each curve refers to the temperature in Centigrade degrees. 


Mean data based:on 67 determinations are given in Table II; the average 
deviation of individual results from the mean was 1.1%. As shown by Fig. 4, 
the data are represented to 1% by the relation 


k X 10° = 29.7 — 0.0687 + 0.11C, 


where TJ and C denote respectively the temperature, °C., and percentage 
Napalm concentration by weight. Changes in concentration, despite their 
marked effect on some of the properties of the gels, produced relatively slight 
changes in the thermal conductivity. 


TABLE II 


THE THERMAL CONDUCTIVITY OF NAPALM-GASOLINE GELS 








k X 10° (cal. sec.-! cm.~! per degree C.) for stated Napalm 








Temp., concentration 
°, 
0% 2.5% 5.0% 7.5% 10% 
50.0 27.0 26.4 26.7 2053 are 
25.0 27.8 28.6 28.4 28.8 29.4 
0.0 29.4 30.0 30.2 30.4 30.8 
—25.0 oe. 31.8 S2:2 32.3 32.4 


—50.0 32.8 33.5 33.7 33.8 34.2 
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A Note on the Temperature Coefficient of Resistance of Unaged 
Tungsten Filaments 


The use of I. C. T. data for the temperature coefficient of aged tungsten 
filaments (7, Vol. 6, p. 136) led to high values for the thermal conductivity 
of known liquids. Determined values for the coefficient of the unaged fila- 
ments used were about 20% smaller than the I. C. T. figure. Two methods 
of measurement were employed. (1) Determinations of the resistance change 
AR caused by a temperature change AT ranging from 2° to 8° C. were made 
with a thinly lacquered filament supported in the constant temperature bath. 
The lack of dependence of the resistance on filament current indicated that 
the filament was closely at the bath temperature. The temperature was 
measured to 0.01° C. with the bath thermometer, or at —25° and —50° C. 
with a calibrated copper—constantan thermopile (three couples). The coeffi- 
cient appropriate to TJ») and Rp» was calculated from ap = AR/AT Ro. 
The resistance measurements were made with the bridge previously described. 
(2) Ro values obtained from the intercepts of the P vs. p/ curves of con- 
ductivity data at bath temperatures from —50° to 50° C. were plotted against 
the corresponding TJ) values. The coefficient at each working temperature was 
calculated from the slope of the curve and the Ro value at that temperature. 

Asummary of the data is given in Table III. The ‘Method 1’ values are the 
means of from three to six determinations. The ‘Method 2’ values are based 


TABLE III 


THE TEMPERATURE COEFFICIENT OF RESISTANCE FOR UNAGED TUNGSTEN FILAMENTS 


Dia | Method a X 104, per degree C. 


50.00 | 62408 


25.00 | 38. 0. 
0.: 


0.00 


—25.00 





—50.00 1 6 + | 
2 0° £004 
| | 


on all the conductivity data, which involved results for at least two different 
tubes at each working temperature. The data are described to better than 
1% by the relation 


a) X 10° = 427 — 1.897) + 0.008875 — 0.0000217%. 
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CHOOSING AN EMULSION AND PROCESSING TECHNIQUE 
FOR DAYLIGHT AERIAL PHOTOGRAPHY! 


By L. E. HowLett 


Abstract 


Experiments are reported that show the relations existing between gamma, 
useful exposure range for aerial photography, maximum resolving power, and 
emulsion speed as determined by the resolving power criterion. Aero Super XX 
and Panatomic X, processed in various ways, were used for the experiments. 
The results are used to deduce the proper choice of emulsion and processing 
technique for ordinary daylight aerial photography where the main requirement 
is the acquisition of information. 


Purpose of the Experiment 


In a previous communication (1) a resolving power criterion was proposed 
for selecting the photographic exposure that, for aerial photography, can be 
considered inversely proportional to the emulsion speed. Such a method 
gives appropriate emphasis to the importance of obtaining information. This 
is generally the principal user requirement in aerial photographic operations. 
[t was shown that, although numerical determinations of speed by the proposed 
method have to be made by a subjective estimate of the limit of resolution, 
good and certainly quite adequate agreement can be attained at different 
times by a number of different observers of varying experience in such work. 

The present discussion reports investigations to determine the relation 
between speed, as determined by the resolving power method, and certain other 
physical characteristics of two emulsions, when several different processing 
techniques are employed. Two emulsions, Aero Super XX and Panatomic X, 
and three developers were studied, D-19b, DK-20, and DK-600. The 
processing time was varied to secure a wide range of gammas. 


Experimental Procedure 


Sensitometric printings were made from a resolving power step tablet. 
On each density step were three identical resolving power targets. A target 
consisted of a series of three line units. The series in respect of size had as a 
common ratio ¥/2. Ina unit the lines were equal in width and to the space 
separating them. The line length was nine times the separation. The lines 
were lighter than the density step forming the background. The logarithm 
of the brightness ratio was 0.2. The exposure time was one-thirtieth of a 
second. The illumination was equivalent to mean noon sunlight modified 
by a Wratten No. 12 filter. The usual controls necessary for accurate 
sensitometry were employed. 

1 Manuscript received November 19, 1947. 


Contribution from the Division of Physics and Electrical Engineering, National Research 
Laboratories, Ottawa, Canada. Issued as N.R.C. No. 1712. 
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After processing the printings of the resolving power step tablets, the maxi- 
mum resolving power attained on each step was read with the aid of a binocular 
microscope. The magnification and illumination were adjusted so that 
reading conditions led to the highest evaluation of the resolving power. 
From the results a curve was plotted to relate resolving power with the 
logarithm of the exposure. 

The shortest exposure, which gave 90% of the maximum resolving power 
attained on the negative material for each particular processing technique, 
was taken as inversely proportional to the speed of the material for that 
technique. 

Discussion of Results 
A very large number of resolving power — log exposure curves were obtained, 


but their complete presentation here would hardly be justified. Accordingly 
a typical one is shown in Fig. 1. The normal characteristic curve has been 
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Fic. 1. A typical resolving power —log exposure curve (A) with the density —log exposure 
curve (B). 


included. Ordinates refer to the step number of the smallest target resolved 
and to the density of the step. The exposure axis is the negative logarithm 
of the exposure in meter-candle seconds. It will be noted that in the sample 


given there is a relatively short range of exposures over which 90% maximum 
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resolving power can be obtained. This seems to be a general condition except 
for the finer grained emulsions at low gammas. The range is the important 
one for aerial photography and, to utilize it fully, accurate exposures must be 
made. 

Repetition of the speed determinations for a specific emulsion and processing 
technique, at widely different times and by different observers, showed the 
same reproducibility of results that was reported previously (1). It is inter- 
esting to record that, early in the experiments reported here, concern was 
caused by obtaining from a series of tests an average log exposure for 90% 
resolving power on Aero Super XX at a gamma of 1.4 which differed by 0.12 
from the results previously reported (1). On investigation the records revealed 
a difference in the emulsion number for the two experiments. Fortunately, 
some of the material bearing the original emulsion number was still available. 
Speed determinations were repeated, and the results demonstrated that no 
discrepancies existed in measurements of material bearing the same emulsion 
number, but that a speed difference on the resolving power basis really existed 
between the two emulsion batches. 

Graphical presentation of the results is the easiest way to show the general 
conclusions that can be drawn from the large mass of experimental data. 

Fig. 2 presents the negative logarithms of the exposures, which correspond 
to 90% maximum resolving power for particular processing techniques, 
plotted against gamma. Speed increases with gamma for both emulsions. 
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Fic. 2. The relation between gamma and minimum log exposure for 90% maximum 
resolving power. 

Fig. 3 shows the maximum resolving powers attained by different processing 
techniques plotted against gamma. _ It will be noted that for each emulsion 
the resolving power is only moderately affected by processing. The extreme 
range of resolving power for Aero Super XX is just about two steps and for 
Panatomic X a little less. The highest values are obtained only by great 





HOWLETT: PHOTOGRAPHIC EMULSIONS AND PROCESSING TECHNIQUES 63 


sacrifice of speed. There is in general a difference of about one step between 
the resolving powers of the two emulsions for the same developer and the same 
gamma. 
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Fic. 3. The relation between gamma and maximum resolving power. 


Fig. 4 gives the range of log exposure over which it is possible to obtain at 
least 90% of the maximum resolving power of each emulsion for each processing 
condition plotted against gamma. With both emulsions the range decreases 
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Fic. 4. Relation between gamma and A log exposure over which 90% maximum resolving 
power can be obtained. 


continuously with increasing gamma for each developer. For a particular 
gamma the range is not greatly affected by a change of developer for Aero 
Super XX, but with Panatomic X some considerable differences are apparent. 
At low gammas the latter emulsion has very long ranges over which 90% 
maximum resolving power may be obtained. This characteristic no doubt 
accounts for its popularity in miniature photography. For the same purpose 
with this emulsion DK-20 will be the preferred developer, since it gives a 
slightly higher resolving power and a higher emulsion speed than the other 
two. 
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Choice of an Emulsion and Processing Technique for Aerial 
Photography 


In choosing an emulsion and processing technique for aerial photography a 
suitable compromise must be made between speed, resolving power, and the 
range of exposure over which 90% of the maximum resolving power of the 
emulsion can be attained. The desirable condition would be to have the three 
factors at a maximum simultaneously, since this condition leads to the greatest 
recovery of information and tolerates greater errors in exposing. Unfortun- 
ately a compromise must be made. Provided resolving power is not unduly 
sacrificed, and providing the range of exposure over which 90% of maximum 
resolving power can be obtained is not reduced below the average brightness 
range of the subject from the air, speed is of great importance. The majority 
of aerial photographic lenses presently available perform better when stopped 
down to f/11 or f/16. High performance aircraft are extensively used in 
aerial photographic operations and short exposures are essential to avoid 
excess image movement unless some form of compensation is provided. 
Consequently, any slight increase in the absolute resolving power of the 
material that is gained by a significant sacrifice of emulsion speed can be 
quite illusory. Indeed, in certain circumstances it might actually lead to a 
smaller recovery of information. For high altitude photography (20,000 ft.) 
it can probably be accepted that the average log brightness range corresponds 
to about 0.75. Consequently, this can be taken as the minimum exposure 
range over which it is essential to obtain 90% of the maximum resolving 
power. At lower altitudes longer ranges are required. 

For general aerial photographic operations the greater speed of Aero 
Super XX makes it preferable to Panatomic X. The maximum resolving 
power of Panatomic X is slightly greater, but it is doubtful whether advantage 
can be taken of the fact in practice on account of the much lower speed. 
D-196 and DK-600 lead to the highest speed with Aero Super XX, and this 
considerable increase of speed does not entail too costly a sacrifice of resolving 
power. Since D-19d leads to a slightly higher speed than DK-608, the prefer- 
ence will be given to it since the other characteristics are, within experimental 
error, the same. No advantage is to be gained from the slightly greater 
speeds obtainable at gammas above 1.4, since the range over which 90% 
resolving power can be attained has been reduced below the acceptable limit 
of 0.75. It will be shown in a later communication that if the information 
is to be recovered by way of paper prints there are other considerations which 
suggest that gamma should not exceed about 1.4. 

Of the emulsions and processing techniques tested, it is clear that Aero 
Super XX developed in D-19d to a gamma of about 1.4 is the most desirable 
combination for aerial photography. Other combinations equivalent to or 
better than this one are probably available. Further studies are being made 
to investigate this situation. 


Reference 
1. Howvett, L. E. Can. J. Research, A 24:1. 1946. 





A SIMPLIFIED BROADSIDE DIPOLE ARRAY! 
By-). H. Brus? 


Abstract 


A broadside dipole array, or ‘billboard’, with a greatly simplified feed and 
matching system, is described. Radiation patterns compare favorably with 
those of the conventional billboard. 


Introduction 


The two most difficult problems inherent in the design of a billboard array 
are: 
(1) To obtain the proper phasing and feed of the dipoles for patterns with 
minimum side lobes and maximum gain; 


(2) To match the antenna, particularly when there are many dipoles placed 
in parallel. 


In the conventional type of billboard array each dipole is fed separately 
from a complex feeder network behind the reflecting screen. The purpose of 
this paper is to describe a billboard array in which this feeder system is almost 
completely eliminated, and the matching problem is simplified, without deteri- 
oration in the radiation patterns. 


The Problem 


It was required to construct an antenna to operate at 90 Mc. per sec. The 
antenna was to be mounted on a tower 200 ft. high, and was to rotate contin- 
uously at 4r.p.m. The array was to have maximum gain and sharpness of 
radiation pattern commensurate with its physical size. Operational conditions 
as applied to the function of the antenna limited the vertical beam width to 
a minimum of 30°. Consideration of the mechanical conditions limited the 
size of the antenna to approximately 14 ft. high and 44 ft. in width. 


The physical dimensions quoted above allowed the possibility of a bill- 
board array of 21 -half-wave (A/2) elements, i.e., seven elements wide by three 
elements high, and horizontally polarized. 


The Array 


To gain simplicity in feed and matching, a vertical stack of three seven- 
element horizontally polarized collinear antennas was used (See Fig. 1). The 
phasing of the individual dipoles in the collinears was accomplished by means 
of quarter-wave stubs, which projected back through the reflecting screen. 


1 Manuscript received June 3, 1947. 
Contribution from the Radio and Electrical Engineering Division, National Research 
Laboratories, Ottawa, Canada. Issued as N.R.C. No. 1707. 
2 Present address: Department of Mathematics, Michigan State College, East Lansing, 
Mich. 
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The spacing between the collinears was one-half wave length. Each collinear 
was fed at its midpoint, and the main point of the feed was at the center of 
the middle collinear. 


PARASITIC REFLECTOR PHASING STUB 


T 


<|~ 


<0 
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Fic. 1. The schematic diagram of the collinear antennas. The screen is placed at a 
distance of \/8 behind the radiators, and the phasing stubs project backwards through the 
screen. 


\ 
| FEED POINTS 


A parasitic reflector was placed at a distance of 25 cm. behind the center 
element of each collinear. The reflector caused the center element to be 
driven harder. This fact, accompanied by the natural taper in the power 
distribution along the array due to radiation, minimized the side lobes. The 
screen, made of horizontal lengths of No. 12 gauge copperweld wire, with a 
vertical spacing of 3 in., was constructed of such a size as to give an overlap 
of X/8 all around the array. 

The radiation pattern (field strength) is given in Fig. 2, The beam width 
at half field strength is 18°, the side lobes are less than 14% of maximum, and 
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Fic. 2. The horizontal radiation pattern of the array. 
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the back radiation is less than 4%. As the array was to be used for common 
transmitting and receiving, the actual radar pattern is effectively the power 
pattern (see Fig. 3). 
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Fic. 3. The effective radar pattern of the array. 


As the impedance of each collinear was in the neighborhood of 500 ohms, 
the impedance of the array turned out to be 165 ohms resistance, with very 
little reactance. This value made it very easy to match into with either open 
wire or coaxial line. The array was matched into a 400-ohm open wire line 
by means of a quarter-wave transformer. 


The bandwidth of the array is nearly the same as: for an equivalent array 
of conventional design, since the error in the phasing of the dipoles, when 
operating off frequency, is nearly the same in both cases. 

The array itself was constructed of 5/16 in. brass tubing, mounted on a 
frame of 2 by 2 in. dressed oak. The brass work was mounted on Johnston 
type No. 65 insulators. The insulators were kept as far as possible from 
high voltage points. The array was sence ted to powers as high as 300 kw. 
peak without any arc-over. 

The feeding of the array during rotation was accomplished by means of 
inductive coupling rings, patterned after an English design. These consisted 
of two flat, tuned rings separated by a Faraday shield. There were no sliding 
contacts. 

Conclusions 


The main features of this array are: 
(1) The simplicity of the feed system. There is only one point of feeding, 
eliminating the necessity of having feeders strung around behind the array. 
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(2) By means of the natural power taper and use of the central parasitic 
reflectors, the spurious radiation is kept toa minimum. The pattern obtained 
is equal in all respects to that obtained using the conventional billboard. 

(3) Using three arrays, of fairly high impedance, in parallel simplifies the 
matching problem. There is only one point to be matched. At that point 
the impedance is such as to make matching simple. Furthermore, as elements 
are added to each collinear, the impedance increases. Thus it is possible to 
build an array whose impedance actually increases as dipoles are added. 





A POSITIVE ION SOURCE! 


9 


By A. J. BAYLy? AND A. G. WARD? 


Abstract 


Positive ions are formed in a low pressure electrodeless discharge which is 
excited by the radio-frequency field of a coil coupled toa 15 Mc. per sec. oscillator. 
The ions are extracted from the discharge through a pumping canal (2 mm. in 
diameter X 12 mm. in length) by a simple arrangement of electrodes. After 
further focusing and acceleration to 50 kv., the beam current at the target is 
500 wa. and the focal spot is 6 mm. in diameter. -Magnetic analysis has shown 
that the beam consists of 50% protons when using hydrogen in the discharge 
and 60% deuterons when using deuterium. Gas is supplied to the source at the 
rate of 15 cc. per hr. 


Introduction 


Protons or deuterons, accelerated to high energies, are used as tools for 
nuclear research, and for other purposes. The usefulness of the accelerators 
producing these high energy particles may be measured in terms of the 
maximum energy of the particles, the energy spread of the beam, the maximum 
current of atomic ions, and other criteria of a similar nature. Since the 
performance of an accelerator is closely related to that of its positive ion 
source, a list of the desirable properties of ion sources that are used in such 
accelerators can be compiled. The source should have a high efficiency 
as measured both by the ratio of the number of emergent atomic ions to the 
number of emergent gas molecules and by the ratio of atomic ion current to 
total beam current. The ion beam from the source must be suitably col- 
limated. The source should give an adequate output of ions, various acceler- 
ators requiring ion currents from a few microamperes to many milliamperes. 
Other desirable properties are ruggedness, long life, and low power con- 
sumption. 

Many positive ion sources have been described in the literature, and they 
can be classified in three main groups. The first group, ‘Canal-ray Sources’, 
use a high voltage, cold cathode discharge to produce ions, the useful beam 
emerging from the source through a small hole in the cathode. Papers by 
Bouwers, Heyn, and Kuntke (2) and Craggs (3) contain data on this kind of 
ion source, and many references to earlier work. In the second group, 
‘Capillary Arc Sources’, a low voltage, high current discharge is passed through 
a constriction or capillary in order to produce a region of intense ionization. 
Positive ions are removed from this discharge either by diffusion through a 
small hole or by an electrostatic extracting field. Sources described by Lamar, 
Buechner, and Van de Graaff (6) and Zinn (9) are good examples of this group. 
The third group of sources are characterized by the use of a magnetic field to 
‘prevent the escape of electrons from the discharge, thus allowing the formation 

1 Manuscript received August 28, 1947. 
Contribution from the Nuclear Physics Branch, Chalk River Laboratory, Division of 


Atomic Energy Research of the National Research Council of Canada. Issued as N.R.C. No. 1703. 
This paper was presented at the Quebec Meeting (May 26, 1947) of the Royal Society of Canada. 


2 Physicist, United Kingdom Staff. 
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of an intense plasma of electrons and positive ions. In sources of this type 
described by Finkelstein (4) and Ardenne (1) a reflection principle is also 
used. Electrons from the cathode are constrained by the magnetic field to 
oscillate between two electrodes at cathode potential until their energy is 
spent in excitation and ionization. In this way the effective path length of 
the electrons is increased and the discharge can be operated at a greatly reduced 
pressure. 

The positive ion source described in this paper does not belong in any of 
the three groups discussed above. <An electrodeless discharge is used to 
produce the positive ions, which are then extracted from the discharge through 
a pumping canal by a suitable arrangement of electrodes. 


Preliminary Experiments 


Experiments on the design of a positive ion source were undertaken as part 
of the work involved in the design and construction of a 50 kv. accelerator. 
A simple glass capillary source, shown in Fig. 1, was used in the first experi- 
ments*. The discharge was produced by a radio-frequency oscillator coupled 
to two external electrodes, and an electrostatic field was used to extract 
positive ions from the discharge through a canal of small diameter. Beam 
currents of 200 wa. were obtained, which contained 20% protons. 
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Fic. 1. Capillary type ion source. Radio-frequency voltage was supplied to the external 
electrodes A, A, to excite a discharge. Positive tons were extracted from the discharge through 
the canal, C, by a steady voltage of approximately 2 kv. between electrodes D, E. Legend for 
aluminum should show crosshatching. , 





* Dr. D. Roaf, former member of the United Kingdom staff, Chalk River Laboratory, began 
the design of the accelerator and made preliminary investigations of this capillary source excited 
by radio-frequency. 
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Radio-frequency excitation of the discharge has two apparent advantages. 
The surface enclosing the discharge region can be made of Pyrex glass, and the 
experience of other workers (5, 6) has shown that this should give a high per- 
centage of atomic ions in the discharge. Since no electron-emitting cathode 
is required, the metal exposed to the discharge is only that associated with the 
extracting electrodes, and the emission failures and indeterminate life of 
cathodes are avoided. 


The properties of this capillary type source were not thoroughly explored. 
Although the beam current was limited by local heating in the capillary, a 
better design for the electrodes would probably have improved the per- 
formance. Recurring breakages of the rather fragile source were trouble- 
some, and the alternative design of ion source shown in Fig. 2 was investigated. 
This design was suggested by the source described by Thonemann (8), which 
was published at this stage of the experiments. 


Fic. 2. Radio-frequency type ion source. The copper coil was connected to the oscillator 
at the points A, A. Positive ions were extracted from the electrodeless discharge in the inverted 
Pyrex bottle by a suitable voltage between Electrodes 1 and 2. A further accelerating voltage 
between Electrodes 2 and 3 was used to focus the ion beam. 


In the initial measurements made on this source beam currents of 300 
a. were obtained, and the proton current was 14% of the total current when 
hydrogen was used in the discharge. With a deuterium discharge the deuteron 
current was 29% of the total beam current. By a variation of parameters 
such as the size and shape of the Pyrex flask and the geometry of the electrodes 
the performance was markedly improved. 
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Final Arrangement 


The arrangement of the ion source used at the present time is shown in 
Fig. 3. The flask in which the discharge takes place is a 600 cc. Pyrex chemical 
reagent bottle. The rim of the bottle is ground flat to give uniform contact 
with the metal electrode. The vacuum seal at this point is made with 
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Fic. 3. Improved model of the radio-frequency type ion source. 


Apiezon Q. A bottle of this shape, with the neck fitting closely on the electrode 
system and a rapid increase of diameter at the level of Electrode 1, gave 
optimum results. A strong discharge occurs near the electrode and a minimum 
of metal surface is exposed to the discharge. 


The exciting coil is mounted concentrically with the axis of the bottle. 
It is made of four turns of } in. diameter copper tubing, and the inside diameter 
of the coil is chosen to clear the walls of the bottle by approximately } in. 
The coil is mounted on two porcelain pillars, and radio-frequency choke coils 
have been wound on these pillars. The ends of the choke coils remote from 
the exciting coil are connected to Electrode 2, Fig. 3. The radio-frequency 
power is supplied through two isolating condensers and a simple parallel wire 
transmission line; connections are made to the coil at the terminals A,A, 
Fig. 3. In this way the 15 Mc. per sec. oscillator is isolated from the 50 
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kv. applied to the source. A push-pull type oscillator is used with an input 
power of 400 w., but it is not known how much radio-frequency power is used 
in exciting the electrodeless discharge. The light emitted by the discharge 
is roughly equal to that from a 15 w. electric light, and a small fan with a 1/10 
h.p. motor gives adequate cooling. 

The dimensions of the electrode system are shown in Fig. 4. The electrodes 
are made of mild steel with the exception of the two duralumin parts shown 
in the figure. Duralumin was used for these parts to minimize sputtering 
from positive ion bombardment. With this source no trouble has arisen 
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Fic. 4. Arrangement of the ion source electrode assembly, the main accelerating gaps, 
G, G, and the target, T. 


from this cause since the discharge appears to remove sputtered material from 
the Pyrex surface as quickly as it is deposited. It is therefore possible that 
an electrode system made entirely of mild steel would be satisfactory. In 
this source the limitation on applied voltage between the upper two electrodes 
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is set by sparking, and mild steel electrodes would raise this limiting voltage. 
The electrodes are spaced with glass disks, and Vinyl-seal resin has been 
used to make the vacuum seals. An assembly jig was used to assure reason- 
able alignment of the components. 


Hydrogen or deuterium is supplied to the discharge from a 3 liter reservoir 
at reduced pressure. It passes through a suitable capillary leak, glass stop- 
cock, and enters the discharge through the inlet tube in Electrode 1, Fig. 3. 
Variation of the pressure in the reservoir is used to alter the gas pressure in 
the discharge. When the operating conditions are normal, gas is used at the 
rate of 15 cc. per hr. The pressure in the discharge, calculated from the gas 
flow and dimensions of the pumping canal, is about 10-? mm. 

Positive ions are extracted from the discharge by a voltage applied between 
Electrodes 1 and 2, Fig. 3, Electrode 1 being positive with respect to Electrode 2. 
A further voltage applied between Electrodes 2 and 3 accelerates and focuses 
the ion beam emerging from the small diameter hole in Electrode 2. These 
voltages, designated as extracting and focusing voltages, are supplied by two 
variable 0 to 10 kv., 10 ma., power supplies. The power supplies are fed 
with mains voltage through a one-to-one transformer with insulation for 
50 kv. between the primary and secondary windings. 


Experimental Results 


Data on the operation of the source under normal conditions are given in 
Table I. The color of the discharge provides a good indication of the behavior 
of the source. Normally the discharge, which fills the bottle, is deep red, 


TABLE I 


OPERATING DATA AND PERFORMANCE OF THE RF ION souRCE (FIG. 3) 


15 cc. per hr. 
~10 X 10-* mm. mercury 


Gas consumption 
Pressure in discharge 


Pressure beyond canal 
Power used in oscillator 
Oscillator frequency 
Extracting volts 

Focusing volts 

Current normally obtainable 
Maximum current obtained 
Atomic ion content 


Ion current density at canal 
Ion current density at target 





~5 X 10-* mm. mercury 
450 w. 

15.5 Mc. per sec. 

1.2 kv. 

7.5 kv. for 50 kv. acceleration 
500 pa. 

750 pa. 

51% from Hz discharge 
57% from Dz discharge 
~16 ma. per cm.? 

~0.65 ma. per cm.? 





and the lines of the Balmer series are very prominent when viewed 
through a spectroscope. If a small air leak is present, the discharge is pale 
pink or white. When the discharge is turned on after the apparatus 
has just been pumped down from atmospheric pressure, it requires one or two 
hours’ operation before the discharge regains its red color. Whengthe operating 
pressure is raised to about twice normal the discharge assumes a ring shape; 
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a reduction of pressure to about one-half normal causes the intensity to decrease 
rapidly, and at still lower pressures only a faint blue glow remains. This 
glow disappears completely when the gas supply is turned off, but fluorescence 
of the flask walls can still be observed. Decreasing the radio-frequency power 
produces similar changes in intensity and the discharge tends to flicker. 


Positive ions are extracted from the discharge by a voltage between Elec- 
trodes 1 and 2 (Fig. 3). The arrangement of the electrodes is such that 
a beam of ions is focused through the small canal in the second electrode. 
For a given intensity of discharge there is an optimum extracting voltage 
corresponding to a maximum beam current passing through the canal, and 
this optimum voltage increases as the intensity of the discharge is increased. 
The relations between extracting voltage, beam current, and intensity of the 
discharge can be explained qualitatively in terms of the position of the plasma 
boundary. This boundary acts asa positive ion emitter, with a current density 
per unit area proportional to the intensity of the discharge. The boundary 
adjusts its position with variation of voltage and intensity of discharge so that 
the positive ion current corresponds to space charge limited conditions, 
When no extracting voltage is applied the plasma projects through the aperture 
in Electrode 1. As the voltage is increased the plasma retracts, and even- 
tually is concave towards this electrode. As the plasma boundary moves with 
voltage, the focusing properties of the ion gun alter, and the operating voltage 
is chosen to give best focusing of the ion beam through the canal. Under 
normal operating conditions, this optimum occurs at 1200 v., and a current 
of about 1.5 ma. is taken from the power unit supplying this voltage. 

When the positive ions pass through the canal they are accelerated by a 
further voltage between Electrodes 2 and 3 (Fig. 3). Adjustment of this 
focusing voltage is used to maintain the focus at the target for variations of 
the main accelerating potential between 10 and 50 kv. Although this electro- 
static lens produces a converging effect on the beam, when it is used in com- 
bination with the two main accelerating gaps (G, Fig. 4) an increase of the 
voltage on this lens causes the beam at the target to diverge. This effect 
arises since the increased voltage of the positive ions entering the main 
accelerating gaps leads to an increase in the focal length. The first crossover 
of the beam after leaving the canal occurs at the target. 

The relative positions of the ion source, main acceleration gaps, G, and the 
target, T, are also shown in Fig. 4. The high voltage is divided equally 
between the two gaps. For measurements of beam current the target, 7, 
was replaced by a Faraday cylinder with appropriate bias voltages and limiting 
stops to avoid errors in measurement due to secondary electrons. 

't has been found that the useful beam current increases from about 200 
ma. at 15 kv. to 500 wa. at 50 kv. if the beam size at the target is maintained 
by appropriate variation of the focusing voltage. The increase in this voltage 
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as the over-all voltage is increased draws an increasing number of ions into 
the beam at the expense of ion current striking the interior walls of Electrode 2 
(Fig. 3). 

The percentage composition of the beam has been determined with a mag- 
netic analyzer. This was placed at the position normally occupied by the 
target, and ion currents were measured in a Faraday cylinder after 180° 
deflection in a magnetic field. The resolution was such as to allow complete 
separation of the mass peaks as far as mass 6, and was determined primarily 
by the beam diameter at the entrance to the analyzer (of the order of 6 mm.) 
and the diameter of the aperture in front of the Faraday cylinder (of the order 
of 1cm.). Flat-topped peaks were obtained as the magnetic field was varied, 
corresponding to one component of the analyzed beam sweeping across the 
aperture in front of the Faraday cylinder. The mean radius of curvature 
in the analyzer was 5 cm., and the analyses were normally made with an 
over-all voltage of 15 kv. 

A series of measurements taken with the ion source operating normally are 
given in Table II. In addition to those peaks listed in Table II, very small 
peaks were observed which corresponded in energy to atomic ions produced 
by the breaking up of molecular ions after they had been accelerated to the 


TABLE II (a) 


TYPICAL ANALYSIS OF HYDROGEN BEAM 











Mass 1 | 2 | 3 | Others 
| 91 | 37 | 51 am 


Peak current, pa. | 





Mass 1 content = 51%. 
TABLE II (0) 


ANALYSES OF DEUTERIUM BEAM 


(Currents in microamperes) 











Mass 2 

Mass 1 2 | 3 4 5 6 content, 
| | | % 

er ae pend Pee 

Run 1 | 4 | 64 24 ag hs 10 62 
Run 2 | 4 | 50% 23 S. ) oe 15 52 
Run 3 | = 7. ae 3 oe | « 144 57 
Run 4 | 44 | 66 | on 19 58 





Average Mass 2 content = 57%. 


maximum voltage. The beam current used in taking these measurements 
was about 200 wa. The higher percentage of atomic ions when deuterium 
gas is used, as compared with that when using hydrogen gas, is a real effect. 


BAYLY AND WARD: A POSITIVE ION SOURCE 


It could arise from the difference in the impurities in the gases, since com- 
mercial samples were used. The ‘hydrogen’ was stated to be 99.5% He 
with O2 as the main impurity, while the ‘deuterium’ was 95% Dez with He as 


the main impurity. 

The percentage of atomic ions in the beam increases gradually as the 
pressure is reduced, but the total beam current decreases more rapidly. The 
normal operating pressure has been chosen to give good percentages of atomic 
ions and a reasonable rate of gas consumption, but the performance of the 
source does not depend markedly on the pressure, and variations of 25% are 


easily tolerated. 

Various workers have reported an improvement in proton percentage when 
gas mixtures are used in the discharge (5, 7). We have tried mixtures of 
hydrogen with helium, argon, water vapor, and oxygen. In all cases the 
color of the discharge changed from red, obtained with hydrogen, to blue, 
and the proton current was reduced. 

The characteristics of various types of ion source are listed in Table III. 
Representative sources have been chosen for each of the three types listed in 
the Introduction. From an inspection of the data given in this table it can 


TABLE III 


COMPARISON OF ION SOURCES 








Capillary arc Reflection 
———_-— type Electrode- 
Type of source Canal ray Metal Glass with less 
capil- capil- magnetic | discharge 
lary lary field 


Lamar, Buechner, Medea Bayly 
and (1) and 
Van de Graaff (6) Ward 


Bouwers, Criaek 
Reference Heyn, and | oo 
Kuntke (2) ‘ 


=20 220 


Pressure in discharge, u 20 
385 135 


Power input, w. 2000 
Canal dimensions, 
diameter X length.mm.| 3 X 5 0.8 X¥0/;0.8 X 0 
Maximum ion current, 
ma. 
Percentage atomic ions 30? 
Maximum energy spread, 
ev. 
H2 consumption, cc. per 
min. at 760 mm. 3.5 
Ratio of protons to gas 
molecules, % 


0.6 . od 1.0 1.5 
20 60 20 


35,000? ~13 ?(Low) 


0.17 








0.07 . . . 0.85 











be seen that the radio-frequency type of ion source has a performance com- 
parable to that of the best of the other sources. It should be noted that the 
values given for energy spread are maximum possible values. In the cases of 
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the Capillary Source with a metal probe (9) and the Radio-frequency Source 
the great majority of the ions in the beam have a relatively small energy 
spread, probably less than 100 ev. 


Conclusion 


The positive ion source described in this paper has operated satisfactorily 
for several months. The ruggedness, simplicity of operation, and good 
performance render it well suited for many of the uses for which positive ion 
sources are required. Improvements in the performance can probably be 
made by a variation of the electrode geometry and an investigation of the 
changes produced by varying the frequency of the oscillator used to excite 
the discharge. 
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GAMMA RAYS PRODUCED IN THE FISSION OF U?*! 
By B. B. KinsEy?, R. C. HANNA’, AND D. VAN PATTER‘ 


Abstract 
This paper describes measurements which have been made on the coincidence 
counting rate between the fissions produced in an ionization chamber and the 
discharges of a Geiger counter. By determining the absolute efficiency of the 
counter, the total energy emitted in the form of y-radiation was estimated to 
be about 4.6 + 1.0 Mev. Measurements have been made on the range of the 
electrons projected from an aluminum radiator placed between the ionization 
chamber and two thin-walled counters. From this experiment, the average 
energy of the radiation appears to be about 2.5 Mev. The result is confirmed 
by absorption measurements in lead. The absorption curve indicates a fairly 
homogeneous radiation and gives some evidence for the existence of a softer 

component with an energy of about 500 kev. 


Introduction 

It was thought that if y-rays are emitted in the process of fission, a study 
of these radiations might yield useful information concerning the mechanics 
of the process. As a first attempt in this direction, the total energy emitted 
in the form of y-radiation, and the hardness of this radiation, has been 
measured. The coincidence method has been used throughout in order to 
distinguish y-rays produced in this way from y-radiation resulting from the 
neutron source and from neutron capture in the sample of uranium and in 
surrounding objects. 


The Coincidence Method 


Fissions are produced by exposing a thin layer of uranium oxide to thermal 
neutrons. The layer is deposited on the collector electrode of an ionization 
chamber connected to a linear amplifier. A Geiger counter is placed nearby; 
the pulses created by the fissions and by the Geiger counter are fed into a 
mixing unit. The coincidences recorded are of two kinds, genuine coincidences 
produced by radiations emitted simultaneously with the fission process and 
chance coincidences. If gs is the rate at which the pulses produced by the 
fissions operate the circuits of the mixing unit, and N is the number of fissions 
produced in the uranium per second, then gy = Ne; , where €; is the efficiency 
of the ionization chamber and the circuits associated with it. If €, is the 
efficiency of production of discharges in the Geiger counter by radiations 
emitted by the fission process, the genuine coincidence counting rate is g = 
Ne;éy. Von Droste and others (3, 7) have shown that the efficiency of a 
Geiger counter for y-radiation is approximately proportional to the energy of 
the radiation, provided that the walls of the counter are thick enough and made 
of material of medium atomic weight. Therefore the probability that a 

1 Manuscript received December 8, 1947. 
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Research Council of Canada. Issued as N.R.C. No. 1723. These measurements were made at 
the Montreal Laboratory during 1944-45. 
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y-ray of energy E,, emitted in the fission process, produces a pulse in the 
mixer is RE, , where k is a constant, determined by calibration, that depends 
on the nature of the counter and on the geometry of the experimental arrange- 
ment. The probability that a pulse is produced by the gamma-rays of all 
energies associated with a single fission is: 


2n,kE, = k2n,E,, 
where 7, is the average number of quanta per fission of energy E,. Thus 
é, = kin,E, = kEy, 


where Eo is the total energy emitted in y-rays per fission. The coincidence 
counting sate is therefore: 

Ge = QskEo. 
Hence 

Eo = 


al 
Ss |S 
~ ° 


The total energy of the radiation may thus be measured in terms of the ratio 
of the genuine coincidence rate to the fission rate when the absolute efficiency 
of the counter is known. The use of the coincidence method enables the 
y-radiation associated with the fission process to be measured separately 
from that resulting from other processes. 


The energy of the radiations has been studied by an adaptation of the 
method of Mitchell and Langer (6). Two counters, with thin walls of dura- 
lumin, were placed side by side and close together with their axes parallel and 
coplanar with that of the ionization chamber. A thick sheet of aluminum 
between the ionization chamber and the first counter serves as a radiator of 
secondary electrons. The triple coincidence rate between the ionization 
chamber and the two counters is measured and plotted as a function of the 
thickness of aluminum absorbers placed between the latter. The thickness 
of aluminum for which the genuine coincidence counting rate is reduced to 
one-half can be used as a measure of the energy of the radiation if this radiation 
is assumed to be homogeneous in energy. 

To determine the degree of homogeneity of this radiation a study has been 
made of its absorption in lead. This has been done by placing lead sheets 
between the chamber and a thick-walled counter adjacent to it. 


Apparatus 


The first experiments were made with a neutron source consisting of a 
mixture of beryllium powder with several curies of polonium, equivalent in 
neutron production to about half a gram of a radium-—beryllium mixture. 
We are indebted to the Health Division of this Laboratory for the loan of 
this source. In subsequent experiments a 1.2 gm. radium—beryllium source 
was used. With the polonium source, roughly equal numbers of neutrons and 
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‘-rays are emitted and it was considered that, with this source, the back- 
ground of y-radiation would be less intense than with a radium-beryllium 
source. However, it was found later that a higher neutron flux could be 
obtained with the radium source for the same background of ‘-radiation 
if the counters were screened from the source by thicker plates of bismuth. 







The neutron source and screening arrangements are shown in Fig. 1. 





NEVTRON SOVRCE 













Fic. 1. Neutron source and screens. 








A large cylindrical cavity in paraffin wax—10 in. in diameter—was formed by a 
paraffin cylinder and cylindrical paraffin plates at top and bottom. The 
neutron source was placed in the center of massive cylinders of metallic 
bismuth. Four inches of bismuth were used to screen the source from the 
enclosure in which the ionization chamber and counters were installed. 
A bismuth cylinder with walls 2 in. thick screened the counters from the capture 
radiation produced by the neutrons in the paraffin and another bismuth disk, 
2 in. thick, closed the top of the cavity. With the radium—beryllium source 
mentioned above, the flux of thermal neutrons in the enclosure was measured 
with the aid of a boron trifluoride counter and found to be of the order of 
1000 thermal neutrons per sq. cm. per sec. 












The ionization chamber used in these experiments was constructed with 
platinum electrodes, mounted in a lead glass tube. A total of 25 mgm. of 
a sample of uranium oxide, enriched in U** by a factor of about 20 to 1, was 
deposited—as uniformly as possible—on a thin platinum sheet with an area 
of about 25 sq.cm. This sheet was rolled to form a cylinder, 8 cm. in length, 
to make one electrode of the ionization chamber. The other electrode con- 
sisted of a cylinder of platinum somewhat longer than the first, and concentric 
with it. After some preliminary heating in vacuo, the chamber was filled with 
pure argon at atmospheric pressure and sealed off. Under these conditions, 
the time taken to create the pulses produced by the fissions is determined by 
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the mobility of free electrons moving in the electric field between the electrodes. 
With 1000 v. between the electrodes, the collection time of the electrons is less 
than a microsecond. To preserve the sharp front of these pulses as they 
are formed at the collector of the ionization chamher, the latter was-connected 
to a linear amplifier of the ‘fast’ type with an upper frequency limit of about 
2 Mc. 

Measurement of the Total Energy 


For these measurements the ionization chamber was supported centrally 
in the enclosure shown in Fig. 1 and contained in an aluminum tube joined 
at its ends to a thin, perforated steel sheet. Four similar aluminum tubes, 
parallel to the first, were disposed at equidistant positions about it, and served 
as supports to locate the Geiger counters. Several copper-walled Geiger 
counters were made up for this purpose, each about 1 in. in diameter with an 
active portion about 8 cm. in length. They were filled with a mixture of 
argon (9 cm. pressure) and alcohol (1 cm. pressure). Owing to the very 
high counting rates which were unavoidable in these measurements, it was 
essential to refill the counters after two or three days of continuous use. , 

The distance between the axes of the tubes containing the ionization 
chamber and the Geiger counters was made as small as possible—4.4 cm. 
The location of the counters and the ionization chamber with respect to each 
other was a matter of considerable importance in securing reproducibility in 
the measurements. To minimize errors which might result from inaccuracy 
in location, two Geiger counters were therefore used successively and the 





results were averaged. 

The measurements fall into two parts: (a) measurement of the ratio of the 
coincidence rate between the ionization chamber and a Geiger counter to the 
fission rate in the chamber, and (0) the absolute calibration of the counter in 


the same conditions. 


(a) THE FissIon-y COINCIDENCE RATE 


In order to verify that the results obtained were independent of instrumental 
conditions, a number of measurements were made with different coincidence 
resolving times and with different potentials across the ionization chamber. 
The counting rates in the Geiger counters were considerable—about 200 per 
sec.—so that it was necessary to correct for chance coincidences. This cor- 
rection was measured directly in each case, by determining the coincidence 
rate between the fissions and the pulses produced by another Geiger counter 
mounted outside the thermal neutron enclosure, screened from it and operated 
by a separate source of y-rays at a similar counting rate. The coincidence- 
resolving time was calculated from the usual formula, g., = 2g;qit, where 
gi is the counting rate in this external counter and 7 is the resolving time. 

With the polonium source, the fission counting rate was nearly constant 


throughout the period of these experiments—about six weeks—at about 1.9 
fissions per second. With the radium source, slight variations in the screening 


‘ally 
ined 
bes, 
‘ved 
iger 
lan 
e of 
ery 
was 


tion 
cm. 
ach 
y in 
acy 
the 


the 
the 
- in 


KINSEY ET AL.: GAMMA RAYS PRODUCED IN THE FISSION OF U*% 83 


arrangements were made from time to time, with the result that varying fission 
rates were used, up to 4.4 fissions per second. Fig. 2 is a typical bias curve— 
one of many—obtained with a discriminator connected to the output of the 


FISSIONS PER M/N-.+/28 


B/AS VOLTS 


Fic. 2. Bias curve for fission pulses. 


amplifier. In order to operate the multivibrator circuits of the mixer unit, 
the magnitude of the pulses must exceed a certain threshold. The gain of 
the amplifier connected to the ionization chamber was adjusted so that the 
natural a@-particles produced pulses lying below this threshold and only 
fissions would work the mixer. The threshold was equivalent to pulses of 
about 10 v. peak, the gain was fixed so that the knee of the bias curve came 
just below this amount. Above it, the bias curve is nearly flat. Thus, the 
efficiency of counting fissions is nearly unity and independent of the gain. 

The results are collected in Table I. The Series A refer to one counter, the 
Series B to another. The errors given are the statistical errors of the counts. 
It will be seen that there is no significant change in the coincidence rate with 
different potentials across the ionization chamber; nor is there any significant 
variation of the coincidence rate with the resolving time. The error in €y 
will depend on the statistical accuracy with which the fission counting rates 
have been measured and on the accuracy of location of the counters with 
respect to the ionization chamber. The average value of €, for all results 
is 3.0 + 0.2 X 10%. It is probably accurate to 10%. 
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TABLE I 
FISSION-’Y COINCIDENCES 


(All counting rates in units of 128) 








Corrected Genuine 


Total coinc. a : 
Counter rate /hr. comer rate/hr. per min. €, X 10 


Polonium—beryllium source 


uae .03 
.30 03 


.24 + .03) 


.29 od 
8 


31 
25 + .03) 0.08 
20 + .02] 1.5 | 0.05 





“ao ‘ 0.66 | 0.03 + 
0.21 01 0.66 | 0.024+ 
0.34 : 7.0 O17 
0.206+ . 0.7 0.033 





sjco ooo soe 
woloco © oo & & 








Weighted mean 





Radium-beryllium source 


0. 22 
0.356 
0.45 





0.27 0. 


Mean 


(6) CALIBRATION OF THE GEIGER COUNTERS 


To a sufficient accuracy, the efficiency of counting of these Geiger counters 
is proportional to the energy of the radiation. Provided that the quantum 
energy emitted in fission is not too different from the energy of the y-radiation 
used to calibrate the counters, calibration at one energy is sufficient. For 
this purpose the 2.62 Mev. radiation emitted in the transition ThC’’—> 
ThD was chosen and the absolute strength of the source determined by 
measuring coincidences between the various products involved in the branching 
disintegration, ThC—>ThD. The results were checked by similar measure- 
ments with the long-period cobalt isotope. 
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About 65% of the atoms of ThC disintegrate with B-emission to form ThC’. 
This is a short-lived body with a period of the order of 10-7 sec. (4). The 
remainder (35%) emit a-particles of lower energy giving ThC”’, a body which 
may be isolated from the thorium active deposit by the method of recoil. 
The calibration was made in two ways, firstly, with the aid of a source of ThB 
in transient equilibrium, and secondly, with a source of pure ThC”’. 


3 X 1077 sec. 
ThC’ 
ee 


ae / 
ThB———ThC 


10.6 hr. 
\ 


3.2 min. 


poet. 
rhe” 


In the first calibrations with subsidiary apparatus, the coincidence counting 
rates between a@- and B-particles has been measured, correcting for the fact that 
not all the B-particles detected follow the transition ThC —> ThC’ —> ThD. 
This measurement gives the rate of disintegration of ThC’ and, hence, 
that of ThC”. The active deposit was then placed in the counter assembly 
in the position normally occupied by the uranium layer and the counting 
rate for y-rays determined. Following von Droste, 70.5% of this counting 
rate was assumed to be due to the 2.62 Mev. radiation. 


In the second calibration, the source of ThC” was placed in the location of 
the uranium layer and its absolute activity determined by measuring coin- 
cidences between a y-counter, used in the fission measurements, and a thin- 
walled B-counter nearby. It was assumed that each disintegration of ThC”’ 
leads to the emission of two y-rays in cascade, one of 2.62 Mev. and the other 
of 0.58 Mev. It is by no means certain that all disintegrations of ThC” 
follow this course, but it does seem established that very nearly one quantum 
of the 2.62 Mev. radiation follows every transition (1). The contribution 
of the lower energy radiations being in proportion to their energies and to 
their intensities, the error involved in the result will be negligible if the 
assumption is made that 0.58 Mev. is the mean energy of these radiations. 


(1) The a-8 Coincidences 

. A source of ThB was deposited on a thin foil of aluminum (about 5 mgm. 
per sq. cm.) by the method of recoil from a radiothorium preparation. A 
small fraction of this source was selected—about a millimeter or two in 
diameter—and secured to a thin aluminum base which was mounted on a 
wire frame placed centrally in the aluminum tube supporting the ionization 
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chamber. For coincidence measurements the wire frame was mounted 
equidistant between the openings of two thin-walled counters each closed by 
a mica window of about 4.5 mgm. per sq. cm. thickness. Of these counters, 
the one was operated as a normal Geiger counter with an argon-alcohol 
filling, the other, facing the active surface, was filled with pure argon and 
used as a proportional counter. The stopping power of the materials between 
the active deposit and the sensitive portion of this counter was sufficient to 
stop all the 4.6 cm. @-particles due to the disintegration of ThC and to allow 
only those to be detected which were due to the disintegration of ThC’ (8.6 
cm. range). The total stopping power between the active surface and the 
B-counter was about 15 mgm. per sq. cm. The £-particles recorded are 
divisible into a hard and a soft variety, the hard type resulting from the disinte- 
gration of ThC” and ThC; the soft type arise from the disintegration 
ThB—>ThC and the internal conversion electrons due to the 0.24 Mev. 
y-ray originating in this transition. The number of these electrons is approxi- 
mately 0.25 electron per disintegration of ThB. A considerable fraction of the 
softer B-radiation will be absorbed in the small amount of material between 
the active surface and the 6-counter. Absorption curves were made with 
aluminum absorbers to determine the fraction of the 8-counting rate which was 
due to the higher energy §-particles. Fig. 3 is a typical absorption curve. 
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50 100 Mgm. jem 















Fic. 3. Absorption of B-particles. 





The first portion of the curve, where the absorption is most marked, is assumed 
to represent the §-particles of ThB and the internal conversion electrons 
referred to above; the flatter, linear portion which follows is due to the higher 
energy electrons, and if this curve is extrapolated back to zero absorption an 
estimate can be made of the fraction of the observed B-radiation which is due 
to ThC and ThC”. From the curve, and to an accuracy of 10%, we find that 
this quantity is 0.65. 
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If, then, NV is the number of atoms of ThB disintegrating per second, NV 
will also represent the rate of disintegration of ThC, because the time allowed 
for the preparation of the active deposit was always sufficient to obtain 
equilibrium between the two elements. If, , 9g, 9,3 represent the counting 
rates of the a-particles, 8-particles, and of coincidences, €a, €s the efficiencies* 
of the counters, and yu that fraction of the B-particles observed which are due 
to the disintegration of ThC and ThC” (wu = 0.65), then:— 


i N €a x 0.65 
Ne€g 


Qug = Ne€a€s X 0.65, 
whence 
N = 49,9/Qap 


(the factor 0.65 representing the proportion of ThC atoms disintegrating to 
form ThC’). The number of disintegrations of ThC” occurring per second 
is given by N” = 0.35N. If Q, is the counting rate in the y-counters when 
the source is placed in the central position normally occupied by the ionization 
chamber, the efficiency for the detection of the 2.62 Mev. radiation is obtained 
from the relation: 

€, = 0.700,/N", 


using von Droste’s result already referred to. 


Some care is necessary in inserting numerical quantities into these equations. 
In the first place, the source is decaying continuously, so that counting rates 
have to be referred to a standard time, which in this case was taken to be the 
mean time of the coincidence run. The observed value of gg is to be corrected 
for the y-counting rate and the background rates due to other y-sources in 
the vicinity and the result referred to this standard time. The coincidence 
rate is to be corrected for the chance rate. The chance rate was measured 
by determining the coincidence rate between the a-particles and the pulses 
produced by an external Geiger counter, operated at approximately the same 
counting rate. The result gives a measure of the coincidence-resolving time; 
it is then corrected to allow for differences between the counting rate of the 
external counter and the @-counting rate in the coincidence measurements 
and for any differences between the a@-counting rates. The y-counting rate 
of the thick-walled counters, Q,, was determined for each of two counters, 
in four distinct positions of the foil holding the active deposit, corresponding 
to the plane of the foil lying in the plane containing its center and the axis of 
a Geiger counter, or perpendicular to it. The mean of these eight measure- 
ments, corrected to the standard time, is taken as Q,. Table II gives the 


* €g refers to the electrons produced in both branches of the transition ThC —> ThD, and ts 
omen be the same for both. 



































88 CANADIAN JOURNAL OF RESEARCH. VOL. 26, SEC. A. 


TABLE II 


a-B COINCIDENCES 
(Counting rates in units of 128 per min.) 





T a8 Chance rate dap N” X : 
Expt sec. Wa 9B dep | - 10-4 Oy Y 
xX 106 Corr. Ext. Corr. Corr sec.~1 a 
A 1.1 120 143 134 0.494 0.078 0.077 0.417 1.88 21.6 i ay 
B 7.0 81 102 91 0.474 0.206 0.196 0.278 1.45 14.4 4.3 
Cc 4.4 49 58 53 0.493 0.078 0.056 0.437 0.32 Sia a4 


results of three measurements. The first two measurements (A and B) 
were made on the same source at times 5.5 hr. apart, and with different 
resolving times. It will be seen that the values for €, agree closely. The 
third measurement (C) was made with a different source and with different 
efficiencies for collection of a- and B-particles. The result for €, is appreciably 
higher. In all these measurements the same value of u is assumed, viz., 0.65. 
The fifth column in the table gives the 6-counting rate corrected for the 
y-counting rate. The correction was made by interposing a thick plate of 
aluminum between the source and the 6-counter. The correction is relatively 
high owing to the high y-ray background in the laboratory. The seventh and 
eighth columns give, respectively, the chance rate—as measured directly 
with the external counter—and as corrected to allow for the different con- 
ditions in the coincidence arrangement. This latter result is subtracted from 
the results of the sixth column to give the figures in the ninth column, which 
are used to calculate N’’. The mean of these three measurements gives: 


€2.62 = 1.9 X 10-?. 
(2) The B-y Coincidences 

The source of ThC” was deposited on a thin layer of aluminum about 1 sq. 
cm. in area, mounted on the wire frame used in the previous measurements. 
The experiments were performed with the frame situated in the middle of 
the aluminum tube used to support the ionization chamber. A hole in the 
wall of this tube allowed §-particles to emerge and operate a thin-walled 
8-counter. As observed above, this method allows the efficiency of the 
Geiger counters used in the fission-7 coincidences to be measured directly by 
using these counters to detect the y-radiation. The source of ThC” was 
prepared by recoil from a source of ThB in equilibrium with its products. 
The purity of the source was tested by plotting decay curves of the activity 
produced; no impurity was detected. 

The procedure followed was to measure the total number of 6 and counts 
and coincidence counts in a period comprising several lifetimes of the radio- 
active deposit (usually 15 to 18 min.). A rough value of the efficiency of 
detecting the y-rays is then given by the ratio of the total coincidence counts 
to the total B-counts. To correct for coincidences due to cosmic rays, and 
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for the background counting rate in the B-counter, the total number of counts 
is again measured in a period of time equal to that of the first period and 
immediately following it. Let pg, and pg, represent the total number of 
coincidences and the total number of 8-particles respectively, recorded in the 
first period, p’,, and p’, the corresponding quantities for the second period, 
If €, is the total efficiency for the detection of the y-rays in this transition, 
then €, is given by the equation :-— 


€é,(1 + 7N) = (Pg, — P'p,)/(Pp — P’g)- (1 + €,/€), 


where 7 is the resolving time of the mixer, €, and €, are the B and y¥ efficiencies 
of the B-counter, and N is the number of atoms of ThC” disintegrating per 
unit time at the commencement of the first period. If the activity has 
decreased to xN at the end of this period, the equation is accurate to small 
quantities of the order of (N r)?, xN €,79 and xN7 , where 79 is the dead time 
of the B-counter*. Three different resolving times were selected and four 
measurements were made in each case corresponding to the four positions of 
the source with respect to the Geiger counter. The results are given in 
Table III; the measured values of €, are corrected to allow for the effect of 
the 0.58 Mev. radiation. In making this correction, it is assumed that the 
counting rates due to the two y-rays are in direct proportion to their energies; 
to allow for this, the value of the efficiency calculated from the figures in the 
fourth column of Table III are multiplied by the factor 0.82 to give the effi- 
ciency for the 2.62 Mev. radiation. The ratio €,/ €, is approximately 7%. 


The errors given in the fourth and fifth columns are the statistical errors 
in counting. The results show considerable deviations among themselves. 
To a certain extent this is due to the low statistical accuracy of the coincidence 
measurements, the total number of coincidences recorded in a run varying 
between 30 and 800, depending on the strength of the source and the geomet- 
rical arrangement. The differences between the mean values of €, for 
different resolving times do not appear to be significant. It is possible that 


* This equation is calculated on the following basis: 

It is assumed that the B-counting rate is big enough to warrant correction for the dead time 
and that the y-counting rate is low enough to ignore tt. The coincidence rate is to be corrected 
for the chance coincidences; the cosmic ray coincidence rates are eliminated from the result by taking 
the differences between the quantities pg, and p'g, assuming that this rate is not seriously affected 
by the dead time. The difference between pg and p'g eliminates the background counting rate of 
the B-counter. To this approximation, a simple calculation gives: 


tate .«,, I= eae ke 
bp — ?'s oS Negro/2.( +x) * (1 €’y/€p)» 


which reduces to the equation given above, where terms of the order of (Nr)*, xN €gto and xNr are 
neglected. The dead time of the Geiger counters was determined by measuring the deviation from 
linearity of the plot of the logarithmic decay of the ThC" source. It was found to be about 150 psec. 
In these measurements €g was usually of the order of 2%. Thus, N being of the order of 104, 
the quantity N€gro/2 is of the order of a few per cent and has a negligible effect on the result. 
For a period of 15 min., x = 0.025 and terms containing x are negligible. 
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TABLE III 


B-y COINCIDENCES 


























7, N X 10-‘sec.-1 TN, Pay — Pav x 30s €:.62 X 10° 
psec. % Pe. — Pe 

(4.2 3 1.73 + .12 1.46 + .09 
0.5 1°55 3 1.47 + .14 1.24 +11 
“< : 2.5 1 1.75 + .10 1.51 + .08 
| 7.0 4 1.73 + .10 1.45 + .08 

Mean 1.42 
Zz 2.1 6 2.03 + .16 1:67 + .12 
a. 1.6 5 1.77 + .17 1.46 + .13 
| 0.7 2 1.36 + .20 1.18 + .16 
{1.0 3 1.20 + .20 1.03 + .17 

7 Mean 1.34 
. +7 14 2.63 + .11 1.98 + .10 
a 2.4 13 2.38 + .13 1.82 + .10 
4.0 21 1.93 + .16 1.33 + .10 
| 4.4 22 3.06 + 112 2:06 + .08 

Mean 1.80 





Mean 1:5: + 15 < 16" 


some genuine coincidences were lost at the lower resolving time; experiments 
which have been made to test this point indicate that the loss is not serious 
until the resolving time is reduced below 0.5 usec. The mean value is 


2 = 135 + 0.15 xX 10-3. 


The measurements of the efficiency for the 2.62 Mev. radiation by the two 
methods do not agree very well but probably not worse than might be expected 
from the crudity of the methods. However, even this low accuracy is suffi- 
cient for rough measurements of this type. Taking the mean of the two 
methods, we obtain 

€é.e = 1.7 + 0.2 X 10°. 


To estimate the total energy emitted as y-rays, this result must be com- 
pared with that for the efficiency of the fission y-rays. The total energy 
radiated is thus found to be 4.6 Mev. 


(3) Calibration with Co®° 

To check on the figures obtained above, the counters have been calibrated 
by measuring the y—y coincidences produced in them when a sample of 
active cobalt is placed in the location of the fission chamber. Deutsch, 
Elliott, and Roberts (2) have shown that the 8-particles emitted by the long 
period radioactive body (5.3 years) are associated with two y-rays emitted 
in cascade, one with energy 1.10 Mev. and the other with 1.30 Mev. After 
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making the usual corrections for chance coincidences, the y-ray efficiency, 
calculated from the ratio of the y—y coincidence rate to the y-counting rate 
in either counter, will represent the efficiency for the mean energy of this 
radiation, viz., 1.20 Mev. If €, represents this mean efficiency, g. and q, 
the coincidence and -counting rates, then: 


€, = 4/g,+(1 + 97), 


where To is the dead time. 


The results are given in Table IV. A and B refer to two runs with the same 
source under slightly different conditions, C to a run with another weaker 


TABLE IV 
-y COINCIDENCES WITH Co*® 


(Counting rates in units of 128 per min.) 





Mean y-rates Coincidence rates X 10-3 
Background ———____—___—— - — ———————— 
rate Counter Counter 
1 2 





Chance Genuine 


196 0.260 .007 | 0.123 + . 0.1374. 
204 0.281 + .007 | 0.137 - 0.144 + 


104 0.107 -002 | 0.034 ‘ 0.073 
85.9 0.061 -001 | 0.022 4 0.039 


source, D to a third source, weaker still. The coincidence resolving time 
was 0.76 + 0.02 usec. The errors given in Table IV are the statistical errors 
of counting. From this table we find that the mean value of €, is 0.78 + 
0.02 X 10-*. This gives 4.6 Mev. for the total energy of radiation in the 
fission process, in agreement with the previous result. 

In the following pages it will be shown that the hardness of the fission 
radiation is about the same as that of the 2.62 Mev. radiation of ThD. Our 
result for the total energy makes it reasonable to assume that no very high 
energy radiation is emitted and that the counter efficiency is indeed linear in 
the energy range concerned. To test the linearity at low energies the effi- 
ciency for the radiations from radioactive gold has been measured in a similar 
manner. These radiations consist of two y-rays of nearly equal intensity 
and energies, 280 and 440 kev. and a weaker y-ray of 70 kev. (5). The 
efficiency was found to be 3.2 X 10-4. This result is about one-fifth of the 
efficiency of the ThD radiation; the ratio of the sum of the energies of the gold 
‘radiation to 2.62 Mev. is rather more than this. 


It seemed possible that a fraction of the measured efficiency of the fission 
radiations was due to recoiling nuclei in the Geiger counter produced by fast 
secondary neutrons associated with fission. A rough estimate of this effect 
gives a total contribution to the efficiency of about 1.10-5,—a negligible 
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amount. To test this possibility, the Geiger counter was operated in the 
proportional region, just below the threshold*. The gas amplification of 
the counter was certainly sufficient to allow all recoil nuclei to operate the 
mixer, because the total counting rate far exceeded any reasonable value for 
the counting rate of the recoils that might have been anticipated. Two runs 
were taken, each over a period of many hours, with 7 = 1.8 usec. In both 
cases, the coincidence rates observed were only slightly in excess of the chance 
rates to be expected; 50% in one experiment and 20% in the other. The net 
result is a maximum possible contribution to the efficiency due to neutrons 
of about 6 X 10-5. This is negligible in comparison with the observed 
efficiency. 

In view of the fact that the uranium layer in the ionization chamber is 
extended over a length of about 8 cm., and the thorium active deposits used 
in the calibration experiments were effectively point sources, it seemed neces- 
sary to determine whether a correction should be applied to the measured 
efficiencies to allow for this difference. This was done by plotting a curve of 
the y-counting rate in the Geiger counter as a function of the position of 
the thorium source along the axis of the tube which serves as a support for 
the ionization chamber. The correction to be applied was found to be less 
than 5% of the total efficiency, and consequently it has been ignored in the 
preceding calculations. The possibility was also examined that scattering 
from the walls of the bismuth enclosure would increase the measured efficiency 
for detection of the fission radiation, as contrasted with that for the thorium 
radiations, for which the calibration work was performed with the apparatus 
outside the enclosure. No detectable difference in the counting rate of a 
Geiger counter was observed when the apparatus with a thorium source in 
position was placed inside or outside the enclosure. 


The Absorption of Secondary Electrons in Aluminum 


To measure the quantum energy of the radiation, the apparatus indicated 
in Fig.4 wasset up. The ionization chamber, J, used in previous experiments, 
was mounted in a steel box together with two thin-walled duralumin counters, 


P Q 


Fic. 4. Arrangement for measurements of absorption of secondary electrons in aluminum. 


G, with their axes lying ina plane. The chambers and the counters with their 


leads were screened electrically from each other. An aluminum plate, P, 


—} in. thick—serves as the radiator of secondary electrons set free by the 


fission ‘y-radiations produced in the chamber. Thin plates of aluminum, Q, 


* The counter was connected to the mixer through a fast amplifier. 
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were placed between the two counters and were used to study the absorption 
of the secondary electrons. In order to give the optimum counting rates 
the counter farthest from the ionization chamber was made bigger than the 
other (1.5 in. in diameter as compared with 1 in.). The experimental proce- 
dure consisted in measuring the triple coincidence rate between the fission 
chamber and the two counters as a function of the thickness of the aluminum 
plates. 


For this experiment, the 1.2 gm. radium—beryllium source was used. The 
steel box was placed inside the container indicated in Fig. 1. The bias curve 
of the fissions produced in the chamber was similar to curves previously 
obtained; the counting rate for fissions was about 3.4 fissions per sec., the 
single counting rates in the Geiger counters, about 400 per sec. Since the 
efficiency for counting triple coincidences in these circumstances is much 
lower than in the simple arrangement used in previous experiments, the 
coincidence counting rates observed were exceedingly low. Most measure- 
ments involved about 100 hr. of continuous operation. To ensure the maxi- 
mum stability of the electrical apparatus, all the electrical units were fed from 
constant voltage transformers. To make certain that power or tube failures 
over the counting periods would not affect the result, the reading of the meter 
recording triple coincidences was plotted against the time. Failures of the 
equipment could be detected by interruptions in the straight lines obtained 
in this way. The triple mixer used in these measurements was similar, in 
principle, to the double mixer previously used. 


The correction for the chance coincidence rate amounted to about 20% 
of the genuine rate in most cases. The chance rate between random counts 
in the three counters was negligible; the high value for the correction arose 
from the fact that some electrons produced externally, or some electrons 
produced by radioactive atoms in the counter walls, could traverse both 
Geiger counters, causing coincidences between them. With no absorber 
between the counters, this coincidence rate was about 30 per sec. The cor- 
rection was determined by measuring the triple coincidence rate between 
these two counters and another situated outside the bismuth enclosure and 
operated with a separate radioactive source. The equivalent resolving time 
was calculated from the formula: 


dh = 29:9 ps7» 


where gq, is the counting rate in the external counter and gg, is the coincidence 
rate between the thin-walled counters. To verify that the electrical con- 
ditions, in fact, did remain constant, the resolving times were measured for 
each absorber and often several times during the course of a run. This 
quantity varied very little from 0.90 wsec. 
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The results are given in Table V. The ‘chance rate’ was calculated from 
the formula above, when q; is made equal to the fission counting rate. 


TABLE V 
ABSORPTION OF SECONDARY ELECTRONS IN ALUMINUM 


(Coincidence rates in units of 128 per hr.) 


Absorber thickness, mgm. Al 
per sq. cm. 0 50 110 145 190 335 


Total coincidence rate 0.036 0.027 0.026 0.0122 0.0173 0.0103 
Chance rate 0.007 0.006 0.005 0.004 0.0035 0.0024 
Genuine rate 0.029 0.021 0.021 0.0081 0.014 0.0079 


The results of these measurements are plotted in Fig. 5. The statistical 
errors in counting are indicated for each point. For purposes of comparison 
a similar curve has been taken for the absorption of the secondary electrons 
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Fic. 5. Absorption of secondary electrons in aluminum. Full line: absorption curve of 
electrons due to the 2.6 Mev. radiation of ThD. 


ejected by the radiation of ThB in equilibrium with its products. In this 
case, the thorium active deposit was surrounded by a lead cylinder 1 cm. in 
radius. The walls of the counters were each about 30 mgm. per sq. cm. in 
thickness. Under these conditions, the electrons observed are due almost 
entirely to the 2.62 Mev. radiation. The fission radiation would thus appear 
to have a similar hardness. 


Absorption in Lead 


To test the results of the previous section, the absorption of the fission 
radiation in lead was examined. The arrangement is shown diagrammatically 
in Fig. 6. A thick-walled Geiger counter was mounted with its axis about 
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9 cm. from the axis of the ionization chamber and screened from it electrically. 
Flat lead plates about 2 in. wide were introduced between the two. The 
coincidence rate between the fissions and the pulses in the Geiger counter 
were studied as a function of the thickness of the lead. The neutron source 
was the 1.2 gm. radium-beryllium source used in previous measurements. 


lead Plates 


Fic. 6. Arrangement for absorption measurements in lead. 


Owing to the relatively low efficiency of detection of the fission radiation 
in this arrangement, the chance coincidence correction is high. The resolving 
time chosen for these measurements was 0.70 usec.; separate measurements of 
this resolving time were made for each absorber to check possible deficiencies 
in the electrical apparatus. The measured values seldom differed from the 
mean by more than 2 or 3%. All the electrical equipment was operated from 
mains supplied through constant voltage transformers. The fission counting 
rate was about 3.7 fissions per sec. The counting rate in the Geiger counter 
was about 360 per sec. This counter was removed and refilled at intervals. 

The results are given in Table VI. The chance coincidence rates were 


measured, as before, by determining the rate of coincidence between the 
fissions and the pulses produced by an external Geiger counter. 


TABLE_VI 


tet § 


ABSORPTION OF FISSION RADIATION IN LEAD 


' Coincidence rates per hour in units of 128 
Absorber thickness, P 


ym. per sq. cm. pe 
8 " wih Genuine 





0.136 
0.098 
0.122 
0.080 
0.071 
0.056 
0.042 
0.043 


eceococeceo 


The results in Table VI are plotted in Fig. 7, together with an absorption 
curve for the radiation of the thorium active deposit for comparison. The 
statistical errors of the points in the fission radiation curve are indicated by 
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the vertical lines. The contribucion of the 2.62 Mev. radiation of ThD 
amounts only to 54% of the total at zero absorption; von Droste gives 70% for 
this figure. Presumably the difference is due to the greater likelihood of the 
softer radiations being scattered into the counter from those parts of the lead 
absorbers situated outside the direct beam of radiation. To test whether the 
bismuth enclosure had any effect on the absorption curve of the fission radia- 
tion, two curves were made with the thorium radiation, one with the sup- 
porting box outside the enclosure, and one with it inside. After correcting 
for the decay of the source, to an experimental error of 1 or 2%, no difference 
could be detected between the two curves at any point. 
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Fic.7. Absorption of fission radiation in lead. 


It will be seen that the average absorption of the fission radiation in lead 
is very similar to that of the hard component of the thorium radiation— 
about 0.043 sq. cm. per gm.—and that this absorption is very uniform. It 
would appear, from this, that the fission radiation is fairly homogenous. The 
previous measurements on the absorption in aluminum of the secondary 
electrons seem to eliminate the possibility that this radiation is, in reality, a 
much harder radiation with a similar lead absorption coefficient. 

The absorption curve of Fig. 7 indicates the presence of a soft component in 
the fission radiation. An attempt was made to separate this soft component 
from the harder radiation by using a thin-walled counter as detector, with a 
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thin layer of gold (5 mgm. per sq. cm.) deposited on its inner surface. The 
gold layer would give a greater sensitivity for the lower energy radiation, the 
thin walls decreasing the sensitivity for detection of the harder component 
by reducing the number of secondary electrons ejected from the walls. These 
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Fic. 8. Absorption of fission radiation in lead. 


TABLE VII 


ABSORPTION OF FISSION RADIATION IN LEAD 


Thickness of absorber, Genuine coincidence rates X 108/128 
in gm./sq. cm. per min. 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 


5.14 
.196 
. 386 
.631 
.772 
085 
.940 
68 

81 
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efforts were unsuccessful; for considerations of sensitivity, it was not possible 
to place the ion chamber and Geiger counter sufficiently far apart to reduce 
the relative number of Compton recoils produced at the walls of the ion 
chamber and passing into the counter. Following these attempts, an effort 
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was made to plot the absorption curve more carefully for small thicknesses of 
lead. The gold-lined counter was used (1.25 in. in diameter), placed with its 
axis 5.4 cm. from that of the ion chamber. In these conditions, the chance 
rates were small in comparison with the genuine coincidence rates. The 
results are given in Table VII. The mean values of these measurements 
are plotted in Fig. 8. There would appear to be some soft radiation; the 
absorption coefficient is about 0.17 sq. cm. per gm., corresponding to an 
energy of some 500 kev. If this radiation really represents a homogeneous 
y-ray, then inspection of the curves indicates that there is probably less than 
one quantum emitted per fission. 

The full line connecting the experimental points in Fig. 7 is drawn to pass 
through the point for zero absorption and to follow the curve of Fig. 8 up to 
4gm.persq.cm. Thereafter, the best fit is made to the experimental values. 
It will be seen that the two curves agree quite well; one point at 5.7 gm. per 
sq. cm. lies rather far from the mean. 


Conclusions 


In addition to the possible existence of soft components, these measurements 
indicate that two radiations with an energy of 2.5 Mev. are emitted in each 
fission process. Unfortunately, the sensitivity of the apparatus is not high 
enough to allow of the verification that two such y-rays are emitted per 
fission, by observing y—y coincidences. In view of the many different ways 
in which the fission process can occur, it is perhaps surprising that the radiation 
is so homogeneous. 
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ANGULAR DISTRIBUTION OF NEUTRONS INSIDE A 
SCATTERING AND ABSORBING MEDIUM! 


By P. R. WALLACE? 


Abstract 


A method is developed for determining the moments of the angular distri- 
bution of monoenergetic neutrons diffusing inside a homogeneous medium. It 
is applicable to the ‘asymptotic’ region, i.e., the region in which surface effects 
may be neglected. In particular, a general derivation is given for the current- 
density relation, involving no assumptions about the relative cross sections for 
scattering and capture. The second moment is expressed in terms of deriv- 
atives of the density, and the analogous procedure for higher moments is indi- 
cated. Moments of the angular distribution relative to arbitrary directions in 
space are defined: general formulae and applications to special cases are given for 
the second moment. Comparison may be made with the results of elementary 
diffusion theory, in which the angular distribution is assumed to be essentially 


linear. 
Introduction 


A considerable amount of work has been done recently on the diffusion of 
neutrons in various types of media, in connection with the study of nuclear 
chain reactions. In so far as the diffusion of monoenergetic neutrons is con- 
cerned, the problem may be approached either through the ‘elementary 
diffusion theory’ (1, 4) of Fermi or through the more accurate ‘transport 
theory’ (2, 3). The latter involves the solution of an integro-differential 
equation and can be treated in general only with extreme difficulty. The 


greatest difficulties arise in connection with the regions close to irregularities 
of the medium (discrete sources, boundaries, etc.). 


The present paper represents an effort to find out as much as possible 
about the asymptotic behavior of neutron distributions, that is to say, the 
behavior in regions sufficiently deep inside a homogeneous medium to justify 
the neglect of surface effects at boundaries and the effects of inhomogeneous 
sources. It is concerned with monoenergetic neutrons only. 


The method used is that of ‘artificial sources’. The idea of this method is 
the following. Suppose one has a mass of scattering and absorbing material 
V; of uniform properties inside which one wishes to determine the asymptotic 
neutron distribution. Outside V; are materials with different and possibly 
inhomogeneous properties. Sources are assumed both in and outside V;. An 
attempt might then be made to compensate for the difference in properties in 
these latter regions by introducing appropriate sources (or sinks) of neutrons, 
and then to treat the whole medium as homogeneous. Now the distribution 
of these sources will, of course, depend on y, the neutron density in the regions 
in question. If the density is known, due to a point source in an infinite 
medium whose properties are uniformly those of V; , then the density every- 
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where may be written as an integral over the artificial sources outside V;. 
In this we are led to an integral equation formulation of the problem. The 
kernel of the integral equation is the point source solution referred to. 

While the method has a number of useful applications in the practical 
solution of neutron diffusion problems, we are concerned here primarily with 
the derivation from it of general results about asymptotic distributions. 

General asymptotic relations have been derived from transport theory 
between density and current (zero and first-order moments of the angular 
distribution respectively) and the higher moments. 

In particular, a general proof is given for the relation 


Current « — gradient (density) 


and the general form of the proportionality factor is given: this has been 
proved heretofore only for the limiting case of weak capture. 


The Transport Equation and Moments of the Angular Distribution 


Let us suppose that we have some general distribution of neutron sources 
in scattering and absorbing media of various properties. We shall designate 
the source distribution as 


q(to, Qo); (2.1) 


this quantity is the number of neutrons being supplied in unit time per unit 
volume and per unit solid angle at ro radiating in the direction of the unit vector 
Q). It is referred to as the density of sources. 


The various media will be characterized by the constants 
1 = total mean free path, 


a = ratio of total to capture mean free path. 






We shall assume that we have a region of uniform propert’es (i.e., constant 
1, @ and q), in the interior of which an asymptotic distribution is established. 
Let us designate this region V;, and the rest of space V.. We shall choose 
the unit of length equal to the value of Jin V;. , 

We describe the neutron distribution by the function W(r,Q) such that 
¥(r, Q)dVdQ is the number of neutrons in the volume element dV whose 
directions are within the solid angle element dQ. 









The transport equation (equation of continuity) appropriate to the problem 
is then 







i= 






where the summation convention of the tensor calculus is used and ‘,2’ indicates 
differentiation with respect to co-ordinates x;. Q (with components Qi :7 = 
1,2,3) is the unit directional vector and Wo is the neutron density in the usual 
sense, i.e., the integral of over all directions 2. (For convenience, we have 
chosen the unit of time such that the velocity of our neutrons is unity.) 
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The terms in this equation have the following significance: 


The first describes neutron efflux from a given volume and solid angle 
element. (It may be alternatively written div (vy).) The second describes 
the removal of neutrons by scattering or capture. 


The first term on the right represents the density of neutrons scattered 


intodVdQ . iF. i 
l l, 

isotropic scattering throughout for simplicity, though the problem is not 
unduly complicated in the case of anisotropic scattering. Since the last term 
represents the source contribution, the equation is a statement of conservation: 


loss of neutrons from dVdQ on the left, gain on the right. 


where /, is the scattering mean free path. We assume 


(2.2) may be written 


OY; +yYy = =" Yo + q+ a(V.) [(: “ in) ¥ 


as 


1 i ~a(r) a) } ices = a]. (2.3) 


U(r) 
where o(V,) is a function whose value is unity in V, and zero in V;. The 
effect of the constant g outside the square bracket on the right hand side is 
simply to add a term g/a@ to y. For convenience we shall drop this term 
throughout. 

It follows from (2.3) that the density in V; is that from a certain distribution 
of (anisotropic) sources in V,. We can designate these sources (the last 
term, in the square bracket, in (2.3) ) by 


x(r, Qo), (2,4) 


which is the density of sources, per unit volume and per unit solid angle. 


The practical applications of the method to date have been concerned with 
certain special conditions. In the first place, so long as the first and third 
terms in the square bracket on the right hand side are concerned, they describe 
anisotropic sources, i.e., sources which do not radiate uniformly in angle. 
An integral equation with an anisotropic point source distribution as kernel 
results. This is very complicated. But in the special cases in which (a) 
qg is isotropic everywhere and (6) the total mean free path / is a constant, the 
anisotropic sources disappear, and the problem is greatly simplified. In the 
case of plane problems, the condition (b) is not necessary, since a change of 
scale may be used to make the mean free path effectively constant. 


The peculiar circumstances in which the method is most effective are those 
in which the inhomogeneities occur in small regions. Whereas the range of 
- the integral in the integral equation is always infinite in the usual Wiener- 
Hopf formulation, it is then small in our formulation and the problem readily 
lends itself to approximate solution by iteration methods. 


However, none of the above restrictions are necessary for the derivation of the 
general results discussed in this paper. 
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At an element of volume dV(r) let the density per unit solid angle dQ, due to 
a point source at fo radiating in the direction Qo in a homogeneous medium, be 


Q(r a to ’ Q, Qo) . (2.5) 

The density in V; is, then, 
y (r, 2) = [ [xe 2) O(r — to, 2, %) dM aV, (2.6) 

Qo ‘e 


and from this the moments may be written down immediately (in terms of 
yY in the region V,). We are interested, however, in discovering relations 
between asymptotic forms of the lower moments which are independent of 
conditions in V,. 

The zero, first, and second moments are: 


y (r) = : [xen Qo) OQ (r — t, Q) dH dV, (2.7) 
29 J Ve 
o> (r) = [ [ x(v Qo) qe (r — fo, Q) dQ dV, (2.8) 
Wy (rf) = : [xr Q) QO (4 — To, %) dQ dV, (2.9) 
No Ve 
where 


2, = [0,9,...9, (ae, 


the integral being taken over the complete solid angle of directions 2. Let 
us then consider the moments of Q. 


Consider the transport equation (2.2) for a point source in a homogeneous 
medium. Let us define the three-dimensional Fourier transform: 


o(k) = [x e* dV(r), 
for which the inverse transform is 


1 

Species —i kr 7 

aoa | oe dV(k) . 

If we take the three-dimensional Fourier transform from r- into k-space of 
(2.2), we get 


$(k, 2, 2%) = ——! eae do (R, %) + &(2, a), (2.10) 


g(r) = 


1—-#.2 4r 
where $(R, Q, Qo) is the Fourier transform of Q(r, Q, Qo), and go is its zero 


moment with respect to Q. 6,({2, Qo) is the angular delta-function with the 
properties: 


62(Q2, Qo) = 0 for Q x Qo, 
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and 


f62(Q, Q)dQ = 4. 


It follows from (2.10) on integration over Q, that 


1 


1 
bo (h, So) = 1 — ik. tank\’ 
1 — (1 — a) ; ) 


and, on multiplying first by Q‘, that 
a, tank 
: _ i-a k ik’ Q;, : 
iy (& Se) = Tee tanh Bt Tak.’ 


(2.12) 


where k = |k|: these are immediate. We wish also to calculate $% . 
Now it may be shown on expressing the components of in terms of polar 
co-ordinates with the vector R as axis, that 


Qi QI 1 [ tank 1 tan—'k a 
i-ei’* an {3 k - a(t ie ) Joi 


tank * 1 “| kik (2.13) 
zoe ke fj’ : ; 


so that 


ij ib ts 
a 


| i fa’ (: : ey 
k 2 # k 


an7l i pd 
we sae | ey} + (2.14) 


2 # k? 


By taking inverse transforms of (2.11), (2.12), and (2.14) we may deduce the 
zero, first, and second moments of Q. Letting u = R/k, we have 


G= (?, Qo) oe am | $0 (k, Q) e7ikur k? dk dQ, 9 (2.15) 


and the higher moments may be written down similarly. 
It should now be noted that, if we are interested only in asymptotic behavior 
in V;, we may ignore the contributions from the /ast terms in (2.12) and (2.14), 
since these give only contributions from neutrons scattered in V, but not yet 
having suffered a collision in V;. Such contributions are part of the ‘surface 
effects’ which we are neglecting. In the other parts, 2 occurs only through 
the factor 1 


1— tk. %° 
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Let us now introduce 


t [mr = + 1 — ik. Q% X(To, Qo) ef*Fo dQdV = f(R) . (2.16) 


We may then, using (2.15) and the similar equations for Q and Q®@ in (2.7) 
to (2.9), write the terms giving rise to asymptotic parts of the first moments 


of y as: 














~ 1 ‘ 1 . 
yo (1r) = 22m) - (R) eee se e~ikur R2 dk dQ, . (2.17) 
" 1 — (1 — a) 
k 
Similarly 
; tank 
7 (1) a(1 — a) k ia rs —iku.r Q 2 
yy; (rf) = ~2am)* f (R) rorreerrer nena inn ku; dk dQ, , (2.18) 
1-— (1 — a) —— 
k 
and 
—a e—iku.r fiftan-k 1 tan~'k 
Vi (9) = samp [FH maa a(1-“y~) ]é 


1-—(1-—a) E 


34 tank 1 tank } 
saat eae inte | Seeeinaas chen ee oe 5 2 AQ C 
+ 5 re (1 ; ) a. | mits | k? dk dQ, . (2.19) 





Some Properties of Tensor Spherical Harmonics 


It is convenient to introduce tensorial generalizations of ordinary spherical 
harmonics as follows: Let us define a radius vector r, = r,v associated with 
a unit vector vy. Then we may write 


" : (- 1)” 1 
T,,....%,) _ n! ry : i 
, Uf Ny ----t, Yy —_ 1 


the subscripts after the comma indicating differentiation with respect to the 
components of fr, . 

We call this the spherical harmonic tensor of order n formed from the 
vector V. 

The components of T7;,....:;,(v) are linear combinations of the ‘surface 

a cos : 
harmonics’ P™ (cos @) ee mg, 6 and ¢ being the polar angles of v. In 
particular, if we take the polar axis as the x-axis, 
Tiyn = P,(cos 4). 


We may then state the following theorem, analogous to a familiar one in 
the theory of ordinary spherical harmonics: 


eer = LaiX(2m + 1)jn(E\CaT +++ i,(U)Tiy--+in(2) » (3.1) 
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(n !)?2” 
(2n) /° 


jn() = Ve Inu (§) ’ (3.3) 


where J,4,, designates the Bessel function of order (m + 4). 


C, = (3.2) 






jn(&) is defined by 








To prove the theorem it is necessary to show that, for any two unit vectors 
u and v, 






CiT 4.0.4, (#) Ty..:-4, @) = Py (8.9). * (3.4) 






Expanding T;,....;, (wu), by direct differentiation, we get 










ty 
Ti,.-, (4) T;,...;, (v) = e \ 1.3, 5----(2n — 1) wi,---u;,, — 1.3.5..... (2n — 3) 


> ii, Mi,---- Mi, + 1.3.5.... (20 — 5) 


) as n 
2 6; ‘ 6; Uj... Uin + ores : or (=) | 
1°2 a6 6S } nN i ly rt, etn ie = 1 ’ 


where the sums are taken over all possible combinations of subscripts in the 
Kronecker deltas. But now all except the first term goes out, all others 
involving Laplacians of 1/r,. Let us now choose the 1-axis in the direction 













of u. Then 


nn ee _, _ (2m)! - a fs _ (2n)! 
Pkcset (wu) i ees (vy) = 2"(n!)3 (— 1) [sox (=) ] ia = 2n(mppat v), 


from which (3,4) follows. 







Calculation of the Density y in the Asymptotic Region 







shall expand f (R) in spherical harmonics: 






1 4 ace / 
dz i (2n + 1) Caf (R) Tiny (4) ; (4.1) 






f(k) = 









where the subscript ‘;,)’ indicates m indices 1;----i,... Then 


fim (k) = [i (R) Tiny (u) dQy. F (4.2) 





Letting £ = 1%, —itu.e me 3(— §)"(2m-+ 1) jalht)CoTm(u)Teto). (43) 












* The equation (3.1) reduces, with the help of (3.4), to the form given in Stratton, Electro- 
magnetic Theory, p. 409 (Equation 56) (McGraw-Hill, New York, 1941). 





1 See Appendix for proof. 
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Putting the expressions (4.2) and (4.3) into (2.17), and using the ortho- 
gonality of the 7's of different order, and the fact that 








. : : _(2m)! 4r 
7 (m v) Tm i) T; “eeee ; 1Q = ee eee ee pF seeee i 
10) f Tons (0) Ti (W) A = ET ig OE (A) 
we get 
1 e 1 
YO). = 5 oou | ————aa 
3 GaP = ei an k 


=.(—1) "(2n+ 1) Crfin (Rk) Tin (2 In (kr). Bdk= 55 =e (- 1)"(2n+1) Ci, Tn) (v) 


” foxy (R)jn(Rr) 
em, kdk. (4.5) 
Ri thiecu tank 





k 


Now from (2.16) it may be seen that f(R) has no singularities for | k| < 1. 
Therefore the asymptotic contribution from the integral in (4.5) will arise 
from the poles atk = + iv,(v <1). Let us write 


In(x) = 5 5 (i (x) + he'(x)I * (4.6) 


ne) = (EHO .¢ 
WE (a) = | Her 6 


hy (x) ° si in | 


where 


(Hankel functions) 


9 . . . . 1 
ho (x) = e-* [ conjugate complex polynomial in ~| . 


{ Jt will be sufficient to prove this result for the particular co- ordinate system in which the 
xi-axts is in the direction of v. Then the only non-zero component of T;,-++ im (v) ts that for which 
ji = ---- = jm = 1, for which the right hand side becomes 


(2m)! 


(m!)? 2™ ian = 1 Pm (1), 





where, in fact, Pn(1) = 
But, by (3.4) and the expression (3.2) for Cn, the left hand side ts 


(2m)! , 
(mi)? 2 P(t. v) T;,---- in (u) dQ. 
Since all the Tj, --- im (u) except that for which j, = ---- = jm = 1 contain the azimuthal angle, 
integration over azimuth gives zero except in this case. In the case j, = --- = jm = 1, however, 


the integral becomes 
41 
7242 = —— 
[? - 2m + 1° 
The validity of (4.4) is therefore established. 


* For Bessel and Hankel functions and their interrelations, see for instance Stratton, Electro- 
magnetic Theory (McGraw-Hill, New York, 1941), pp. 358, 359. 
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The last integral in (4.5) may be broken up into 


= \A » (2) 
7 fin (RAS oo — k°dk he ; Fin (RAS Oo kdb. 
1-(-a— i~<twaes 


-2 —o 


hb 


We deform the first of these integrals about the upper half-plane, and the 
second about the lower. 

Since the integrand has branch points at —z and i the contours will be 
semicircular with cuts down the imaginary axis to these points. The contri- 
butions from the semicircles will vanish as the radius becomes infinite. The 
integrals along the edges of the cuts will vanish at least as e~" as r —> ©, 
and will therefore be negligible for large r compared with the contributions 
from the poles at —iv, iv,(v <1). Taking only the asymptotic (pole) 
contributions, we get 


xi [=e fin (iv)h (ivr) 


vr—aea 


= rj [POT feu (— tvher(— bvr) 


e's 


where tv is the zero of E — (1 — a) = | in the upper half-plane, i.e., 


tanh~p . 1 Ss (4.8) 


= eo ’ 


-v 1-—a 
where /, is the scattering mean free path. 
Therefore finally the asymptotic density is 


Yo (as) = cai 1 me) DC — 1)"(20 + 1)CaT in (v) 


vy 
[fim (tv) hn? (ivr) + fim (— tv)hy?(— ivr). (4.9) 


In the particular case of a point source, 


f = f/4e = fan (4.10) 


only the first term in the series (4.9) remains, and we get 


1 v(1— v*) 
OD GRD genet icc ecccentetaieane en eee 
¥ 2mrr (v2? — a) (1 — a) = (4.11) 


Calculation of the Current yj” in the Asymptotic Region 


Let us start with (2.18): 
tan—'k 
ey ee . . 
= am JF @- ae ee 
1 — (1 — a) E 
Cf. Bothe (2, 3). 
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We may write this 








tank 
+) (1 — a) os k ik Q 
=_ — 2: dlldtnccvigenie: sii tmaesipiaeaabetinaiaitiaatiieaciaastinaalicialaicimtaie —tRUT 
i= vi 2(2m)8 f (R) ce e dkd =e (5.1) 
i-(l-—a ; 


where V; indicates the i** component of the gradient operator. The asymp- 
totic part of the bracket may be calculated exactly as was Y ) in the 
preceding section; it is in fact found to be 


= yo (as) z 
p2 


so that 


(1) (@s) __ a {10 
1 = a Ve eh. 











Thus we have the general form of the asymptotic current-density relation. 


a 1 esate : ‘ : 
As v —> 0, a —> which is the value taken in the elementary diffusion 


Z ’ 
theory. Expanding, 


a i 4 108 , z 
S=5[1+3e+ie +--+]. (5.3) 


Calculation of {7 in the Asymptotic Region 
The asymptotic part of Yj” may be calculated from (2.19): 


a _l1-@ gee | - 2° — tank ] . 
4 “~ 2(2ar)8 re cl aces tank|2L k k? k oi 
k 


at ~ 
ne E 1 (1 ae or) _ i tan =] wan RdkdQ,, . 





k 2 =k 


We may write 








1—a e—ku.r sl tank 1 tank 
=e pg ee a Ae ee 46 BALAO 
2(Om)8 a ee oe E i(1 E ) ae | ui; Haka Q, 
1-—a e~*u.r Sit tan—!k 1 tan—'k 
=—V ii ceding ica asc aie Ea sania’ Seige IB wicbe. see, ameter O ‘ 
"7 | 2(27)8 Pe —<—=al 3 (1 k ) 2 sk Jaxaa. | 


(6.1) 


The asymptotic part of the integral in the bracket may now be calculated 
by contour integration to be 


; = (1 aa 3) yy) (as) (6.2) 


1 
v 
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Similarly 
e~iku.r 1 = 1 ( © tank 
tank 2] k k? k 


— a) ; 
a. ae <) YO Cas) | 


—*)| RdkdQ, 


1—(1 
(6.3) 


It follows, then, that 
we. 3a a 
za(t =) vy ' 


asymptotically 
(2) (as) 1-&) yo (as) 6:; 
By 


i 
Let us now define 42(m) as the mean square cosine of the distribution with 
f ’ 
(6.5) 


respect to the direction defined by the unit vector ”, fixed in space 
oe "Fae 


definition 
w(n) 
1 Hye (as) 


Substituting from (6.4) we find that 
2 yO cs (8) On? 


Me(1) 
indicates differentiation in the direction n 


0 
where — 
on 
For very small capture cross section 
3a 4 
=> SS 1 + —= py? 
v 15 


eg 
45 


2,/,(0) (as) 
—. (6.7) 


1 
on? 


and therefore 
2 iS ee 
to" + 15 yo Cs) 


In the first approximation, this is independent of m, the angular distribution 
other hand, for large capture, a, y —> 1 and 
(6.8) 


being isotropic. On .. - 
1 eyo (as) 
win) = yO (ae) ont 
If we wish to evaluate y? relative to a direction which is not fixed in space 
(as, for example, the radial direction in spherical polar co-ordinates), we must 
0? 0? On; O 
= +5 -’14>->.- 6.9 
on? ” Ox, Ox; (6.9) 


put 
nn; 
wT 8x,0X; 








110 CANADIAN JOURNAL OF RESEARCH. VOL. 26, SEC. A. 


oe? 


In spherical polar, co-ordinates, where n; = = this becomes 55. If yy @) 
depends on r alone, 
Tyo (as) — py (as) | (6.10) 
from which it follows that 
1 eyo (as) 2 1 oy (as) oe 
PYO s) gp i= = poy or hws (6.11) 


where yp is the average cosine of the angular distribution relative to the radial 
direction. Therefore 
=o a 3a 1 — 
= %4+(3-1)2a. (6.12) 
Let us consider two special cases here. 
(a) If we have a point source, 
yom) = QO, 


and since 





- 
ll 
Bee 
a 

Is 
“ee” 
=> 
<q 
a 
= 
ll 
|Z 
| 
S12 
$S|+ 
ae 
3 


we have 





1 — 1 1 
man eee = 6.13 
ar B vr vr? ( ) 
Therefore 
SF ce Oe she cw GN ee 
= o+ ( i 1) ( + i) (6.14) 


We note that p? approaches a constant as vr —> ©, i.e., at distances large 
compared with a diffusion length. This may be seen to be an instance of 
a general result. 


(b) If we are in an enclosed spherical region with sources outside, 








yo (as) = C sinh vr : 
r 
Then, as in (a) above, 
aoe ah : coth pr (6.15) 
ar = Tap Or ; , 


which for small pr is 


1 Re 
oi ial (1 “a r) 
Therefore 


eos tee 3a he i eee ee 22 
“a ee 1) 5(1 s™) state 1) 7, Se) 


aaa wl eat > -—= § 
which indicates how the angular distribution departs from isotropy G = 3) 


as we go away from the center. 





WALLACE: ANGULAR DISTRIBUTION OF NEUTRONS 


Evaluation of 2 in General 


We may wish to express y? relative to directions which are not fixed in 
space, and in non-Cartesian co-ordinate systems. For this purpose we must 
express in generalized co-ordinates nini'V;;Y@ and V?y@), The 
tensorial forms of these quantities are n'n’D;D; yO“ and D'D; yO @, 
where D; is the operation of covariant differentiation. 


Now 


niniD,D YOO = nini | cas ‘an {# 5 yO) (as) (7.1) 


Ox'dx! 


{#} being the Christoffel symbol of the second kind.f 


Let us take a co-ordinate system in which x! = s is distance measured normal 
to the surfaces of constant density (i.e., in the direction of the current). 
x?, x’ may be arbitrarily chosen co-ordinates in the surfaces of constant 
density. Let the indices a, 8 take the values 2, 3.* Now y “*) depends 
only on s. From these facts, and the definition of the Christoffel symbol, it 
follows that 

n'niD,D YO =  « 2 - 1 n°n® ag 9 | you ‘ (7.2) 

ds? 2 Os OS 


Consider next 
Vay (as) — at D, Dy (as) — pry (as) ** (7.3) 


Calculating, the Laplacian for the particular co-ordinate system in question, 
one finds 


yp (as) = E de 5a aB rat yi (as) , (7.4) 


Substituting for the second derivative with respect to s in (7.2) above, we find 


_2 alee a on 
ete) = 5(1 7 =) ” sapien (3 1). 


1 op Ogg OY =} 4 1 jaye Map ee. (7.5) 


1\2 2 4f,(0) (as) — sali 
[ one oy 2° Os Os _ Os Os 
But 


_@ ay (as) 


a 


py Os 


= = myo, 
and therefore 
-, = a 
w(n) = 5(1 -4) +5 
1 04.3 — 1 
1)2 gat See 7.6 
[ne{t +55 4 as # 2a om 
+ See, for example, McConnell, Applications of the Absolute Differential Calculus. (Blackie, 
1931) (pp. 140 ff.). 


* The line-element is do® = ajjdxidxi = ds* + aggdx%dx®, 
** McConnell, p. 151. 
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Let us suppose that we may choose the co-ordinates x? and x? in such a 
way that 






(7.7) 


Qep = Sap» 
















where a,, do not depend ons. Then 


ie 1 Be 

ciate B ee ap Seiya OIE Sos ee 11 sea 
~ = aa. = atiq ala = 

2a Os as -< as (a¥ai; un) 


(7.8) | 


as’ 















t int te 2 a x 12 
ae ae ee egg ee [1 (n')*] . (7.9) 


With the aid of these facts we may write 
az _*)\ 4 4. 
w(n) = 5 (1 9) + (3 1) 
[mye{t + 2 cm)} — © 2 — (ony wn). .10) 
\ as as ; 


Now n' is the cosine of the angle which our direction m makes with the current. 
Designating yw for the forward direction by fo, 


w(n) = n'. Mo (7.11) 


wa = 11-2) 448-9 
[cope -— n'| — Hag d= ms | (7.12) 


Further we may write, by virtue of (5.2) 


and 


(7.13) 









If we introduce 
Sp SS Sy - (7.14) 


7H (0-8) +408 — for - Bae 


So 
+ {1 — (”! je log yo]. (7.15) 


An interesting special case is the following: if we calculate y? for a direction 
perpendicular to the current, we get, from (7.15) and (7.11) 


a ad Oe 
B= 5(1 ni) 


which is completely independent of geometry, and of the particular position 
chosen. 
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We might also note the result, that at a sufficiently large distance in units 
of 1/v, i.e., at a sufficiently large number of diffusion lengths, outside a closed 
surface containing all the sources, we get 


a(n) = 5(1 - 4) + 3(¥ a 1) (n'y, (7.16) 


v 
° ° er ‘ — a 
which is again independent of geometry. This reduces to pw? = 7 for n! = 
1, i.e., 2 in the direction of the current. 
Relation Between the Moments 


Multiplying the transport equation (2.2) by {2 and integrating, we get, in 
the asymptotic region, 


yi es — Yi; : (8.1) 
Introducing 

aa (2) ij 
so that eS oe 

wn) = (wn n;, (8.3) 
we get i 

yor = —[(w)*po],, 

= — (MY, — (WH) (8.4) 


ee Qa oe . . . ° . 
lf (y*)4 = 7 6" in the asymptotic region, this would reduce immediately to 
(5.2). However, it may be seen in fact that this is not so: (uw) is not 
diagonalized, or even constant in space, in the asymptotic region. 

To show this we use the asymptotic form of (4) derived from (6.4): 


ies. 1 at t £3 (0) (as), 
(uw)? = 3 (1 — a/v)6% + Pr (33 - 1\ oe (8.5) 


We see that (#°)‘’ varies with position in the asymptotic region and does not 
reduce simply to a diagonal tensor. The second term on the right hand 
side of (8.4) is not negligible. It is only by combining the two terms (as may 
be_immediately verified) that we get the relation (5.2). 


Higher Moments 


We note finally that the methods used in calculating moments up to the 
second can easily be extended to moments of higher order. In calculating 
the Fourier transforms of the higher moments, it is necessary to evaluate 


(see (2.10)), or alternatively, 


/ rr an (9.1) 


“4—&. g 
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This is essentially the problem of determining the tensorial generalizations of 
the Legendre functions of the second kind, just as T‘-----™ is the generalized 
Legendre function of the first kind (Legendre polynomial). 

Designating the m‘* moment, JQ*"Q' YdQ, as y™ir---im, the m™ moment 
relative to the direction of a unit vector n is 


un(n) = ni----m,,Yrirr im, (9.2) 


° e . . ° th ~ ° 
will be expressible in terms of derivatives of Y up to the m™. Conclusions 
th . a ° e . 
about the n"' moments relative to definite directions may be drawn as in the 
cases considered above. 
APPENDIX 


Expansion in Tensorial ‘Spherical Harmonics’ 
The function f(R) = f(ku) may be expanded in the form 
f(ku) = iz S, (2n +1) Sf (eu) Pr(u. u')dQ , (A.1) 


(see Rverly, Fourier’s Series and Spherical Harmonics: (Ginn & Co., 1895) 
p. 211). But using the result (3.4), it follows that ; 


f(k) = i En (2m + 10a [SFR Ten (ud Qu] Tea (1), (A.2) 


from which (4.1), (4.2) follow. 
It should be noted that if we attempt to expand f(k) in terms of the spherical 
harmonic tensors 7;,), the coefficients f(,) will not necessarily be unique. 


(s + 1) @ +2) 


components of 7,,) (after taking account of 





For of the 


2 
symmetry in all pairs of indices), only (2m + 1) are independent, because of 
n(n — 1 : 2 ee 2 
the fact that the ais =) contractions with respect to two indices of T(,) vanish. 


There are an equal number of independent coefficients f(,).. But if the sym- 
metry and contraction conditions are applied also to the f(,»)’s, they also 
become uniquely determined, and have the values (4.2). 
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